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ABSTRACT 


A study on the effects of rotation in the areas of vibrations, 
dynamic stability, optimization, and aeroelasticity is presented. The 
governing equations of motion for the study of vibration and dynamic 
stability of a rapidly rotating deformable body are developed starting 
from the nonlinear theory of elasticity. Through some selected exam- 
ples, some common features such as the limitations of the classical 
theory of elasticity, the choice of axis system, the property of self- 
adjointness, the phenomenon of frequency splitting, shortcomings of 
stability methods as applied to gyroscopic systems, and the effect of 
internal and external damping on stability in gyroscopic systems are 
Identified and discussed. Then the experience gained from this gen- 
eral study is profitably applied to three specific problems . In the 
first problem exact vibration frequencies in two perpendicular planes 
of a rotating uniform beam inclined at an arbitrary angle to the axis 
of rotation are calculated. Also, a stability analysis of this beam 
is performed with the angle of inclination and rotational speed as 
parameters. Existence of several interesting bands of instabilities 
is shown. Optimum depth and thus the optimum mass distribution of a 
rotating beam with a given fundamental flapping frequency and rota- 
tional speed is calculated in the second problem using a numerical 
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iterative method, based on a transition matrix approach, to solve the 
resulting two-point nonlinear boundary-value problem. It is shown that 
even in the presence o£ rotation the optimal structure is characterized 
by the condition that specific Lagrangian density is constant. In the 
last problem the effect of steady coning angle and internal damping at 
the blade flapping hinge on whirl flutter stability in prop-rotors is 
investigated. The developed mathematical equations and their analyses 
are applied to two wind tunnel models studied in previous investigations. 
The results indicate that these new parameters have a marked influence on 
stability, and that they may even change the critical flutter mode. The 
destabilizing effect of damping at the blade flapping hinge is discovered 
to be analogous to the effects of internal damping in gyroscopic systems. 
It is shown that this internal damping might help to explain why some 
previous investigations show good correlation between theory and exper- 
iment and some do not. 
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I . INTRODUCTION 


1.1 BACKGROUND AND OBJECTIVES 

Engineering science provides many examples in vibrations, dynamic 
stability, optimization, and aeroelasticity involving the effects of 
rapid rotation. Some examples are the rotating shaft, the rotating 
disk, the rotating propeller blade, the spinning space vehicle, and 
the spinning planet. These bodies may be of one-dimensional or two- 
dimensional or three-dimensional configurations; may be rigid or de- 
formable; may or may not suffer external kinematic constraints; and • 
may be subjected to both conservative and nonconservative external 
forces . The term rapid rotation means that the additional forces ex- 
perienced by the body due to rotation are of the same order as those 
acting on the body in the absence of rotation. Most of these problems 
are understood in a piecewise manner, and satisfactory analytical 
treatments have been given in the literature. However, several un- 
resolved discrepancies remain in the literature, particularly in the 
areas of vibrations, dynamic stability, and aeroelasticity of de- 
formable bodies. Thus, the subject of effects of rotation covers a 
broad spectrum of both applied and academic objectives, methods and 
applications even when the topic is restricted to structural dynamics. 
Hence, one can recognize multiple stages to a study. The stages 
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considered in this dissertation are: 1) the development of the equa- 

tions of motion of an unrestrained^ deformable, and rapidly rotating 

body, and a systematic reduction of these equations to specific cases; 

\ ■ 

2) a general study of vibration and stability of a rapidly rotating de- 
formable body; 5 ) consideration of some common phenomena such as fre- 
quency splitting, orthogonality between the modes, the effect of choice 
of axis systems, and the effect of internal and external damping on 
stability; and 4) application of the experience gained in stages 1-5 
to three problems of practical engineering interest. 

The first problem considered is a slender beam, capable of bending 
in two planes perpendicular to its neutral axis and simultaneously ro- 
tating about an axis inclined at an arbitrary angle to its neutral axis 
The angle between the neutral axis of the beam and the plane perpendicu 
lar to the axis of rotation is called coning angle, and the beam is re- 
ferred to as a coned beam. The exact vibration frequencies of this 
beam are calculated at various coning angles and rotational speeds. 

The principal motivating factor for this investigation is the fact 
that both Coriolis and centrifugal forces appear due to coning angle.- 
Another is the need for an exact vibration, dynamic response, and 
stability analysis of the rotor blade with coning angle, in view of 
the use of precone in several v/STOL designs. 

The second problem concerns application of control theory methods 
to the mass optimization of structures having dynamic constraints in 
the presence of rotation. The model considered for this study is a 
rotating cantilever beam with given fundamental flapping frequency. 
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The motivating factor for this investigation is a more general one of 
avoiding unnecessary weight, which increases overall weight and in turn 
increases the direct operating costs of v/STOL prop-rotor configurations . 

The third investigation is aimed at better understanding of \^irl 
flutter in prop-rotors. Whirl flutter is an aeroelastic instability 
that can occur on a flexibly mounted aircraft prop-rotor/power-plant 
combination, and involves coupling of inertial, elastic, centrifugal, 
Coriolis, aerodynamic, and damping forces. Whirl flutter is a design 
consideration in every v/STOL aircraft configuration which makes use of 
rotors or propellers for lift or thrust. The primary factor that moti- 
vates this study is the lack of agreement between theory and experiment 
in some specific cases . 

1.2 SCOPE OF THE THESIS 

The common factor in all the problems studied in this thesis is 
rapid rotation. Two types of mathematical models — the elastic con- 
tinuum model and the concentrated mass and stiffness model — are used 
to represent the structure of rapidly rotating deformable bodies. A 
"viscous” representation is adopted for both internal and external 
damping. 

This thesis comprises seven chapters and one appendix. In Chap- 
ter II, using an elastic continuum model and the nonlinear theory of 
elasticity, the equations of motion of a rapidly rotating deformable 
body are developed. Then these equations are linearized about a steady 
state condition in Chapter III, and the equations for free vibration 


- 3 - 



are obtained. Though the explicit need for nonlinear elasticity theory 
is well known in the static theory of elasticity in order to derive the 
equations for deflection of a beam subjected to axial loads, it has ap- 
parently not been used explicitly in deriving the free vibration motion 
equations for rotating beams simply because it is easier to avoid com- 
plexity by treating centrifugal forces by an "effective applied load" 
concept. As a result, the explicit need for nonlinear elasticity 
theory in the presence of rotation is not even mentioned in the liter- 
ature. The effective applied load concept is well suited for rotating 
beams, but for a general elastic body it is not easily implemented. 
Hence, the explicit need for the nonlinear theory of elasticity in the 
presence of rotation is emphasized. With regard to stability, simple 
examples are used to illustrate shortcomings of the kinetic method of 
stability as applied to gyroscopic conservative systems, and the con- 
sequent need for the introduction of imperfections to extract all the 
possible instabilities of the systems. The advantage of using a body 
fixed rotating coordinate system to detect the instabilities due to 
imperfections is emphasized. Also, the effects of internal and ex- 
ternal damping in gyroscopic systems are discussed. Finally, vibration 
and stability of a beam at various rotational speeds are studied with 
coning angle as a parameter, and some interesting instability bands are 
found . 

Chapter IV is devoted to the calculation of exact vibration fre- 
quencies and mode shapes of a rotating beam with coning angle. These 
results are used to validate the application of Galerkin*s method to 
this class of problems . 



Least weight design with dynamic constraints is the problem of 
Chapter V. A rotating beam with given fundamental frequency and rota- 
tional speed is optimized for least weight using optimal control theory 
methods to derive the necessary conditions for an extremal. A nvimeri- 
cal solution technique, called the transition matrix approach, is adapted 
to solve the resulting two-point boundary-value problem. The numerical 
results appear to be somewhat unusual. Optimality of the structure is 
characterized by the constancy of a specific Lagrangian density which 
includes the potential energy density of the equilibrium centrifugal 
forces. This observation is consistent with the corresponding one for 
nonrotating structures . 

A whirl flutter study in prop-rotors in the propeller mode of oper- 
ation, aimed at better correlation between theory and experiment in some 
specific cases, forms the contents of Chapter VI. In this problem the 
power-plant-nacelle-propeller combination is represented by a concen- 
trated mass and stiffness mathematical model. It is shoxm in Chapter 
III that the coning angle of the beam changes the blade flapping fre- 
quency, and that internal damping has a destabilizing effect on sta- 
bility in gyroscopic systems . Given the importance of flapping fre- 
quency and internal damping on whirl flutter, two additional parameters 
— coning angle of the blades and internal damping at the blade flapping 
hinge — are introduced in the formulation. The equations of motion are 
derived by a Lagrangian approach. Quasi-steady blade element theory is 
used to generate aerodynamic forces of the windmilling prop-rotor. A 
linearized version of these equations is applied to two wind-tunnel 
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models from previous investigations, and the whirl flutter stability 
boundaries are obtained. The results indicate that these parameters 
have a marked effect on stability. It is also shown that internal 
damping might even explain why some investigations show good correla- 
tion between theory and experiment and some do not. 

In a concluding chapter are discussed some limitations of the con- 
tinuum mathematical model for analytical treatment of the problems in 
the presence of rotation. Some common features of problems involving 
effects of rotation are summarized. Possible extensions for this in- 
vestigation are suggested. 
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II. EQUATIONS OF MOTION OF A RAPIDLY 


ROTATING DEFORMABLE BODY 


2.1 INTRODUCTION 

The quantitative theoretical treatment of any class of engineering 
problems requires a physical idealization, appropriate laws of physics, 
a resulting mathematical model, mathematical equations, and their anal- 
ysis . In this chapter the governing equations of motion for the study 
of vibration and stability of a three-dimensional, unconstrained, de- 
formable, and rapidly rotating body are derived by considering the body 
as an elastic continuum. The derivation is straightforward, and appears 
implicitly in the literature (Ref. 2.1 for example). It is reviewed in 
this chapter with multiple objectives; 1) to use these general equa- 
tions in the following chapter to extract general inferences concerning 
the influence of rapid rotation on the equations of motion, and on the 
frequencies and mode shapes of a rotating deformable body; 2) to obtain 
the equations of motion of some examples considered in the following 
chapters as a special case of these general equations; 5) to make ex- 
plicit the necessity of employing nonlinear elasticity theory to study 
vibration in the presence of rotation; and h) to maintain consistency 
in notation. 
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2.2 EQUILIBRIUM EQUATIONS 


Let us consider a body of total mass m as shown in Fig, 2.1. 

This body may undergo displacements with respect to an orthogonal body- 
fixed system xyz with the origin at o ^ the center of mass of the 

" ^ ' y — ^ ^ 

body in the undeformed state. We shall denote by - ix + jy + kz 

the position vector relative to axes xyz of an element of mass dm 

in the undeformed state, by lu(x,y,z, t) + jv(x,y, z,t) + iw(x,y, z,t) 

— ^ ^ ^ ^ 

the elastic displacement vector of the element, and by r=r^+q=i| 
+ ^ the position vector relative to axes xyz of the element of 
mass in the deformed state. All these vectors are referred to the same 


set of rectangular Cartesian coordinates xyz . The variables x,y, z,t 
and are called Lagrangian and Eulerian variables, respectively. 

The components of q need not be infinitesimally small. The absolute 
position of dm in the deformed state is given by + r^ + q , vhere 
is the position vector of the center of mass of the body in the un- 
deformed state relative to the inertial axis system x'^y^'z'" 

If the body is acted on by surface tractions per unit of area des- 
ignated by the vector and by body forces per unit of volume desig- 

nated by the vector R , then the equations for linear momentum G 

X 

and angular momentum ^ are given by 



( 2 . 1 ) 


( 2 . 2 ) 
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FIG. 2,1--Rapidly rotating, thre^-dimensional, unconstrained, 
and deformable body - r^ , gpsition vector of mass 
point before deformation; q , elastic displacement; 
Tq , position vector to the center of mass o; Ic' , 
position vector to the mass point after deformation; 

, Steady angular velocity of the body. 



where suffix ^x' Indicates that the corresponding quantities are re- 
ferred to the initial state* and ^ X ( ) » Here^ Tt 

at 0t 

is the angular velocity of the body expressed with respect to the xyz 
system* If we introduce the relation 




= r + r 
o 


( 2 . 3 ) 


into Eqs. (2.1) and (2.2) we obtain 




= mV + 






■ ■ ■ ^ 1' ^ mm ^ 

H = I . n + I'' . u 
o o o 


ii(/' ('u X ^ * fjf (7x 5?) dv. 


( 3 . 5 ), 


Here ^ is the linear velocity of the point o expressed in the xyz 
system^ and and are inertial tensors given by 


I = fff (r^l - r ?) p dV 
O JJJ u U X J 


V 

o 


= Iff 


r , + q I - q q) dV^ (2.6) 


ql-rq-q„ 
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rs/ 

and I is the idem or unit tensor. While obtaining Eqs . (2.4) and 
(2.5)> we made use of the fact that o is the center of mass of the 
body. 

Applying the law of linear momentum (see Ref. 2.2), we get 

-f/f ffK ( 2 - 7 > 

''x 

Similarly, the law of angular momentum leads to 

^ = fff ^ ^x> * If V (2.8) 

V s 

X X 

Equations (2.7) and (2.8) govern the motion of the body, but they pro- 
vide no information concerning the internal or elastic response. The 
equations that govern the elastic response are referred to as elemen- 
tal equations, and they must be obtained by considering the stress and 
strain within the body. 

According to Novozhilov [2.1J, deformable bodies are divided into 
two types: massive bodies and flexible bodies. In the former type, 

the extension is of the same order of magnitude in all three of its 
dimensions; and in the latter type, the extension in one or two direc- 
tions is small compared to that in remaining dimensions. The choice 
of the elemental equations depends on the type of the body. The clas- 
sical equations of elasticity are applicable to all massive bodies and 
to a particular class of flexible bodies (such as rods, plates, and 
shells). In particular, they may be used in the study of bending and 
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torsion of thin rods when they are not subjected to the simultaneous 
action of axial forces, and the angles of rotations of the sections are 
small compared to unity and at the same time small compared to strains. 
Hence^ for rapidly rotating bodies, particularly when flexible, the 
classical theory is no good, since these bodies are subjected to equil- 
ibrium stresses, and the angles of rotations of the sections are not 
small compared to strains . Hien, the next question concerns the type 
of nonlinear theory to be used. This was discussed by Milne [2.3]. 

The present writer agrees with this discussion, since the equations for 
most engineering type problems can be derived with certain assumptions. 
The assumptions that are appropriate for this nonlinear theory are as 
follows : 

(1) No body couples, no couple stresses. 

(2) Elongations and shears small compared to unity. 

( 5 ) Infinitesimal strains. 

(h) Linear stress-strain law (Hooke’s law). 

It is to be noted that there is no restriction on the magnitude of dis- 
placements and rotations. The appropriate elemental equations of equil- 
ibrium, compatibility, boundary conditions, and stress -strain Law can be 
deduced from Ref. 2.1 as 


3 — - - ^ 

dt^ 


[a . (I + 


= 0 


(2.9) 


ij 


O'. . 


E.,, = e. 


j + 2 


( 2 . 10 ) 
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and 


. (I + . n = 

X X X 

(2.11) 

. 2^E + . ? 

(2.12) 


where 



(2.15) 


The suffix 'x* , which indicates that the quantities are referred to 
the initial statej will be omitted here onwards for simplicity, since 
the areas and volumes do not change due to deformation in view of the 
assumption that shears and elongations are small compared to unity. 


2.3 SPECIFICATION OF BODY AXES 


In the absence of kinematic constraints the elastic boundary value 
problem is the Neumann problem. The solution of such a problem is arbi- 
trary to the extent of a rigid-body displacement and rotation. This ar- 
bitrariness can be resolved with the choice of the body axes. A thorough 
discussion on choosing the appropriate body axes was presented by Milne 
[2.U]. In all, there are three possible choices: namely, attached, 

mean, and principal axes . Attached axes are specified by the condition 
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that the rigid-body displacements and rotations are zero. Mean axes 
are chosen in such a way that at every instant the linear and angular 
momenta of the motion with respect to the body axes are identically 
zero. The basic requirement in the case of principal axes is that the 
tensor I' should be diagonal. It was also illustrated by Milne [2.3] 
that the conditions for both mean axes and principal axes will be iden- 
tical for shapes which are typical of aeroplanes, in which transverse 
displacement relative to a plane or line contributes the main deforma- 
tion. We shall use the mean axes in all our discussions. These axes 
are given by the conditions 



It is to be noted that the first of Eqs. ( 2 . 14 ) inq>Ues that the 
origin of the mean axis system xyz is always at the center of mass 
of the body. Substituting Eqs. (2.IU) into Eqs. ( 2 . 7 ) and (2.8), we- 
obtain the overall body equations as 


^ ^ fff ^ 

V S 


(2.15) 



+ I" . r?) 
o ^ 


~ jjy + JTCr X^)ds (2.16) 

V s 



There is no change in the elemental Eqs . {2.S)> (2.10), (2.11), 
and (2.12) of the body. 

We note that the motion of the body is described by a "hybrid" set 
of equations, since Eqs. (2.I5) and (2.l6) are six ordinary differential 
equations, and Eq. (2.9) leads to three partial differential equations. 
Further, all these equations are nonlinear. These equations form the 
bases for mathematical analyses of the vibration, dynamic stability, 
and aeroelastic behavior of an unconstrained, three-dimensional, rapidly 
rotating and deformable body. The equations of motion of some of the 
examples considered in the latter chapters are obtained as special cases 
of these general equations. 
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III. VIBRATION OF A RAPIDLY ROTATING 


DEFORMABLE BODY 


3.1 INTRODUCTION 

In the previous chapter we developed the equations of motion of a 
rapidly rotating, three-dimensional, unconstrained, and deformable body 
for the study of vibration, dynamic stability, and aeroelastic behavior. 
In this chapter we shall deduce the linearized equations of motion of 
vibration of the same body in a vacutma and shall present a brief dis- 
cussion concerning the magnitude and character of the influence of ro- 
tation on the equations of motion of vibration, and on the frequencies 
and mode shapes. Then, these discussions will be followed up by several 
examples, though all of them are not new, in order to extract some com- 
mon features of the effects of rotation. These features include choice 
of axis system, the property of self-adjointness, the phenomenon of 
frequency splitting, shortcomings of stability methods as applied to 
gyroscopic systems, and the role of Internal damping on stability. 
Finally, vibration and stability analysis of a rotating beam with 
coning angle as a parameter will be presented. This analysis is 
believed to be new. 

The experience gained in this chapter, particularly on the effect 
of coning on the frequencies of the rotating beam and on the influence 
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of internal damping on stability in gyroscopic systems, is applied to 
the whirl flutter problem in Chapter VI. 

5.2 EQUATIONS OF VIBRATION IN A VACUUM 

The most general model for the vibration problem is an uncon- 
strained, three-dimensional, rapidly rotating, deformable body, whose 
motion is assumed to be stable. A small perturbation motion about some 
stable steady-state condition can be called vibration. As mentioned 
earlier, the elastic problem is of a Neumann type, and hence the solu- 
tion is arbitrary to the extent of a rigid-body translation and rota- 
tion. This arbitrariness was resolved by properly fixing the reference 
axis system, which is the mean axis system in our case. These axes will 
translate = constant) and rotate = constant) with the body. 

Motions, \diich satisfy the equations of motion when no external forces 
act on the system with respect to the mean axes, might be called gen- 
eral motions. Of all these motions, only those involving steady non- 
rotating translation = constant , if = 0) of the mean axes are 

normally called free vibration motions. All other motions involving 
steady translation and rotation = constant , ?f = constant) of 

the mean axes will be here onwards termed as vibrational motion with 
rotation, and the associated frequencies and modes will be called fre- 
quencies with rotation and modes with rotation. 

In the absence of external forces, the equations of motion for 
vibration of a steadily rotating, uniformly translating, unconstrained, 
and deformable body can be deduced from Eqs . (2.I5), ( 2 . 16 ), (2.9), 
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(2.10), (2,11), and (2.12). They are with respect to the mean axes as 


follows: 


d 

— (m^) = 0 

dt 




P dV 


= 0 




dt 


- ^ . [a . (I + ^^)] = 


cr. (^ + ^^).T? = 0 


(5.1) 

(5.2) 
(5.5) 
(5.i^) 


Equations (2.10) and (2.12) are unaltered. 

The above set of equations are nonlinear. We shall linearize 
them about the steady state condition. Let us use the suffix '1’ 
for the characteristics of the perturbed condition (?f= , 

suffix ’o' for the characteristics of the steady-state condition 
and the same quantities for the perturbations. Then, the disturbed 
motion will be of the form 


qi = + q 





(J + (J 

O 


t. = If + if 

1 o 


(Cont'd) 
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and 


r + q + q 

it * 


(p^)^ = (p>) 


Application o£ Eqs. (5.1), (5.2)^ (5.5), and (5.4) at both steady- 
state condition and the disturbed condition leads to 


X [(if + ??) X ( t ^ + q^o)] 




+ q X (if X ^) [ dV^ J = 0 


div [cr. (l + ^^ + ^?)]+^ . (a 

X^O X ' ' 


d q 


[cr* (l + ^q +^^)j - if+cr = 

o o ^ 


(5.8) 
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and 


^ . (I ^ = p ^ + q„) (J.9) 

dt 

It should be noted that the equation governing the overall linear 
translational motion is automatically satisfied^ whereas that of the 
angular motion couples the rigid body motions and the elastic deforma- 
tions. If = 0 , it can be seen from Eq. ( 5 . 6 ) that the coupling 

between rigid body motions and elastic motions will be eliminated. 

We shall make two more asstimptions : 1) the disturbances are small, 

which enables us to discard the terms a • grad "q from Eqs. (5*7) and 
( 3 . 8 ) and simplifies Eq, (2.12) with the use of Eq. (2.10) as 


cr 


2{xe -h iXv • q 


( 5 . 10 ) 


and 2) the difference between the geometry of the initial undeformed 
state (natural state) and the undisturbed steady state (ft^ = constant , 
^ = constant) is comparatively small; this allows us to omit the terms 
grad from Eqs. (3*7) "" (3.9)* 

After the above simplifications, the linear equations or the var- 
iational equations of the system can be written as 


j2— > 

d q 


dt‘ 


- ^ . a + ^ . ( a . ^q) = 0 


(3.U) 
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and the boundary conditions as 


^ ^ ^ i ^ y V 

cr*n + a^*Vq»n = 0 (5*12) 

Equation (5-6) remains unchanged, and the equilibrium stress tensor 
cr^ can be obtained from Eq. (5*9) 

^ . a = p {?? X X (r^ 4- ^ )]] 

o o ^o 

= p (?f?f - . (r ^ q ) (5.13) 

^00 O U O ^ / 

Substituting Eq. (5*10) into Eqs . (5*H) and (3.12), we obtain 

a q - ^ -4 -4 

p — - + Q x(« xq)+2ft X — 

c O O O 

St dt 

- \i^q- (X + jj) ^ - ^ . ^7) = 0 (3.li^) 

(^ . if . + a .^q.if = 0 (3.15) 

s o 

where suffix represents the symmetric part of the operator ^ , 

Using the principle of virtual displacements, Eqs* (5.11) and 
( 5 . 12 ) can be combined and expressed in integral form as 

j2— > 1 

^ d q 

^.a-*f^.(cr — pr * &? P dV 

° dt^ J ^ 

- J'J’ (^ • n - <y^ ♦^^•nf). dS = 0 (3.16) 

S 
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where 6q^ is the virtual displacement. Equation (3.16) can be 
written in a compact integral form if we take into account that 


( 5 . 17 ) 

V 

^ . 8^) = (^ . a) . 8^ + a : (%-^f 

= (^ . cr) . 8i^ + a : ^^8^ (5. 10) 


T^ere A is any vector and represents transpose. While simpli- 

fying Eq. ( 3 . 18 ), the symmetric property of the stress tensor Z is 
used. Consequently^ we obtain 


///» 


3 q 

— ^ X (if X + 2if X ~ 

[ at^ ° ° ° 0t 



dV 


V 


*//Po ° (3.19) 

V 

To ensure equilibrium, the rotational Eqs . (3.6) must be satis- 
fied. The physical explanation for this is more involved. A discus- 
sion on this was presented by Milne [3.I]. ihe necessary conditions 
for steady rotation of the vibrating body are: 1) mean and principal 

axis must coincide, and 2) the diagonal element of r which corresponds 
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to the axis of rotation must be zero. These two conditions ensure that 
the rotation about the appropriate principal axis of the reference con- 
figuration is stable. For the elastic body^ the principal axes are 
given by the condition that the inertia tensor should be diagonal . 

When the overall motion is stable^ the elemental equations of equil- 
ibrium ( 3 . 1 ^) and the associated boundary conditions given by Eq. (5 *15) 
govern the vibrational motion of the rapidly rotating, unconstrained, 
and deformable body. 

Let us examine Eq. (5. 19) term by term to understand the effect of 
each individual term on the frequency of the body. If ^ = 0 , Eq. 

(5.19) reads 




^ . Sjdv *fffa : 


q dV 


= 0 


(5.20) 


or in the differential form as 


n 2-> 

S q 


— - - (X + p) ^ = 


St 


(5.21) 


and the associated boundary conditions are 


♦ q ) + 2pn . ^q 


(5.22) 


The boundary-value problem represented by Eqs. (5.21) and (5.22) is 
self-adjoint according to Courant and Hilbert [5.2]> since it arises 
from a functional with a homogeneous quadratic integrand. Knowing 


- 25 - 



the problem is self-adjoint^ we can state that the eigenvalues 

2 2 2 

X (X = - CD , CD being the natural frequency) are real, and the 
eigenvectors are orthogonal. It is to be noted that the symbol X 
is used both for the eigenvalue and for one of the material constants . 

A thorough discussion on these eigenfunctions and their use in obtain- 
ing the approximate solutions of dynamic response problems was pre- 
sented by Bisplinghoff and Ashley [ 3 . 3 ], Bolotin [ 5 .^], and Courant 
and Hilbert [3*2] among others. Also, it should be noted that when 
= 0 , the overall Eq. (3*6) of rotational equilibrium is decoupled 

from the elemental equations. 

Now, let us return to the discussion on vibration of the rapidly 

rotating body. The constant rotation adds three additional terms 

to the elemental equations, which can be seen either from Eq, (3.l4) or 

from Eq. (3*19). The centrifugal inertial term [it x x ^1 shows 

o o 

that the natural frequencies (n = 1,2, ...) will be altered. The 
change in the frequencies of a rotating shaft is due to this term. 

Also, this term may give rise to some constant terms, such as the one 
that appears in the equation of motion of a beam rotating about an axis 
inclined at an arbitrary angle to the axis of the beam. The coriolis 
term 2^^ X (Sq7^t) leads to the result that, as shown later, the 
linear eigenvalue problem governed by Eqs . (3.l4) and ( 3 .I 3 ) is not 
self-adjoint. Consequently, the modes will no longer be orthogonal. 

An example of this case occurs in the whirling vibration of a shaft 
with heavy disks attached to it. Another effect of the above two 
terms either individually or collectively is the phenomenon of 
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frequency splitting, which will be discussed later. Another interesting 
aspect of these two terms is that they may cause static instability, as 
in the case of a symmetric shaft, or they may cause self-excited vibra- 
tion as in the case of an unsyometric shaft. Further, the coriolis 
term does have a marked influence on the effect of damping on stability, 
as will be seen in the shaft problems. The stability of some simple 
systems involving these terms will be discussed in later sections. 

The third additional term, ^ • (cr^ • ^7) > changes the fre- 

quencies of vibration. The increase in the frequencies of the heli- 
copter rotor blade with increase in rotational speed is an example of 
this effect. 

When the frequency spectrtim of the nonrotating structure contains 
multiple roots, these roots may be split during rotation into a dis- 
tinct series. For example, a nonrotating elastic shaft has an infinite 
number of natural frequencies in any plane passing through its axis. 

Upon introducing rotation, each natural frequency is split into two; ' 
these frequencies are called whirl frequencies. In another example. 
Backus and Gilbert [5.5»5»6] predicted the frequency splitting for the 
earth’s first two spheroidal modes and verified their prediction, with- 
in experimental accuracy, by seismographic measurements of response to 
the violent Chilean earthquake of May, i 960 . Ihis fact was brought to 
the attention of aeroelasticians by Ashley t5*T]. This phenomenon of 
frequency splitting depends, as shown through examples later, on the 
nature of the system and the coordinate system used. 


- 25 - 



Earlier, we stated that the differential equations and the asso- 
ciated boundary conditions of the vibration problem of a rapidly rota- 
ting body, in general, are not self-adjoint. Now, we shall demonstrate 
this fact. First, we shall assume 

^ ^ , X = + iix> (5.25) 

where ^ is a vector function. Substituting Eq. (5*23) into Eqs. 
(5.1^) and (5«15)> we can write the resulting equations in compact 
form as 

L(J) = x^(y) (3.24) 

\(^) = 0 (n = 1,2, ...,6 ) (3.25) 


where 



) 


(X + p) • a j . ^ 






^ + 


^ X X ^ 2^ X 


-» 


(3.26) 

(5.27) 


Bn(7) 


x^ . 


^ • n + 2pv ^ 


XI + 0“ 




n 


(3.28) 


The boundary-value problem represented by Eqs. (3.24) and (3.25) is 
said to be self-adjoint if, for any two arbitrary comparison vectors 
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^ and ^ , the following identities are satisfied: 

ffft^ . I(?j) dv dV (3.39) 

V 

JJj" ^2 . N(^^) dV = JJJ ♦ N(^2) dV 0 , 30 ) 

V 

Since the tensor differential operator L and the boundary condi- 
tions given by Eq. (5.25) constitute the Euler-Lagrangian equations and 
the natural boundary conditions for the variational problem of finding 
stationary values of a functional with a homogeneous quadratic integrand^ 
it is self-adjoint according to Courant and Hilbert [5.2]. However, the 

/w 

linear tensor operator N contains a nonsyrametric Coriolis term. As a 
result of this, Eq. (3 .50) is not satisfied. Thus, the boundary -value 
problem is not self-adjoint. It is concluded that the vibration modes 
are not orthogonal and cannot exist independently of each other. They 
merely represent the possible harmonic modes of vibration of a rotating 
body. Consequently, it is incorrect to speak of them as normal modes 
of vibration. As stated earlier, they must be referred to as vibration 
modes with rotation. 

A direct analysis of the eigenvalues of a nonself -ad joint boundary- 
value problem is extremely difficult. The possibility of using approxi- 
mate and numerical methods is considered. The most effective approxi- 
mate method is to express the components of the displacement vector as 
a series in a system of coordinate functions (called also as character- 
istic functions or eigenfunctions) which satisfy all the boundary 
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conditions. Application of Galerkin's method, as described “ty Bis- 
plinghoff and Ashley t3.53 several others, reduces the problem to 
one of solving a system of ordinary differential equations in the co- 
efficients of the series . 

Now, there is a physical argument behind the use of normal coor- 
dinates in representing any distortions of the rotating body as the 
sum of distortions in its principal modes of nonrotating body. Another 
issue is the application of Galerkin's method in which convergence 
is not conclusively proved, particularly to nonself-adjolnt boundary- 
value problems . These two issues are more involved, and hence we will 
not go into the details. However, some discussion on these issues can 
be found in Bishop [3*8] Bolotin [3 A], respectively. Our rationale 
to justify this approach is based on a comparison of exact results (if 
not available, experimental results) with those obtained by Galerkin's 
method. A problem, vibration of a rotating beam with a coning angle, 
is solved both by Galerkin' s method and by an exact method (in the 
numerical sense), and the results are compared in the next chapter. 

The solution of Eq. (5.19) will be sought in the form 

00 

7 ^ ^(x,y,z) |^(t) ( 5 . 31 ) 

n=l 

where ^^(t) are the functions of t to be determined, and ^(x,y, z) 
are the characteristic functions of the nonrotating body. Substituting 


- 28 - 



Eq, ( 3 . 31 ) into Eq. (3.19), we obtain the equations of Galerkin' s method: 


M t + 
nn ti 


(tx> 

NR^n o 


M I + 2^ i ^ • 
nti n X ^ r o 


— > 
m 

rn 


r=l 
(Skip r=n) 


- D «/o • 


. Qi 

nr o 


r=l 


where 


E l a : 

r o r n 

r=l 


0 



(n = 


( 5 . 32 ) 


(3 03) 
! 


Comparing the above Eq^ (3*32) with the corresponding one^ Eq* (59) 
of Ref* 3*T> shows that they differ at two places. Firsts is the ad- 
dition of the last term in Eq. (3*52)> which is missing in Eq* (59) 
of Ref. 3 •lx since potential energy due to centrifugal force was not 
considered* The other one is the sign before in Eq. (59)^ which 

should read as positive* Comparing Eq* (3*19) with the corresponding 
one, Eq. (8) of Ref. reveals that the second term of Eq. (8) 
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should be a vector cross product [2if x instead 

of a dot product before the integral sign. Also, Milne replaced the 
gradient operator of the last term of Eq. (5.32) by its antisymmetric 
part, assuming that the effects of rotations of the elements of the 
body are larger than the strains. 

3.3 STABILITY THEORY 

Stability theory is a subject which is of importance in all 
branches of physics. Many questions concerning this theory are still 
unanswered. Some of the controversial issues are: the definition of 

stability, linearization of the equations, correlation between the 
problem to be studied and the methods of stability analysis to be used, 
the role of internal and external damping in gyroscopic systems, the 
use of normal modes of a nonrotating structure for the rotating one, 
and the effect of model truncation. A review of the entire theory is 
out of place in this section, the main purpose of which is to find cor- 
relations between the problems to be studied and the methods to be used 
in order to identify some common features of rapidly rotating, defor- 
mable bodies of engineering interest. However, extensive discussions 
can be found on the concepts of stability in Ziegler [3.9], Leipholz 
[ 5 . 10 ], Stoker [ 5 . 11], Bolotin t3.^J, Hermann [3-12], Bryan f3.13], 
and Melrovitch t5.U] in addition to several others, m this section 
we shall focus our discussion on two main points; 1) choice of axis 

system, and 2) shortcomings of kinetic and static methods of stability 
as applied to gyroscopic systems. 
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Basically the available approaches to dynamic stability analysis 
are knovn as the Liapunov direct method or Liapunov* s second method, the 
dynamic or vibration or kinetic method, and the static methods (imper- 
fection, equilibrium and energy) . The word dynamic stability encom- 
passes all the stability problems analyzed with the aid of Newton's 
equations of motion or by any equivalent method. The idea behind 
Liapunov's approach is to answer the stability question by utilizing 
the differential equations of the system without solving them. The 
method involves devising a suitable scalar function known as the 
Liapunov function in the phase space or motion space and using this 
function along with the motion equations to test the system stability. 
This method is valid for both linear and nonlinear systems . The main 
problem with this method is the construction of the Liapunov function, 
which is not unique, and indeed there is a large degree of flexibility 
in its selection. In the early stages, its applicability was restricted 
to systems governed by ordinary differential equations. In recent years, 
however, this method was extended to elastic and aeroelastic systems' 
governed by partial differential equations by Wang [3-l5]^ and Parks 
[ 5 . 16 ], and to systems governed by hybrid (both ordinary and partial) 
equations by Meirovitch [5,17>5*l83. 

The kinetic instability approach involves an examination of the 
roots of the characteristic equation if the system is governed by 
ordinary differential equations [for instance, Eq. (5.52)]; and an 
examination of the eigenvalues of the boundary -value problem if the 
system is governed by partial differential equations [for example. 
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Eqs. (3.24) and (3.25)1. Also, damping forces, either viscous type 
or structural type, can be included in these equations as shown in 
our following examples. The eigenvalues may be real or complex. 

According to Eq, (3.23), if |a is negative and 0 0 , the 

motion is damped simple harmonic; if X is entirely real and nega- 
tive, the disturbed motion settles back to steady motion. In either 
case it is stable. If p, is positive with oj / 0 or X is entirely 
real and positive, the displacements increase exponentially with time. 
Hence, the motion is unstable. If X = 0 , the displacements linearly 

increase with time in the case of nongyroscopic systems, and a static 
divergence boundary is defined. This is the approach used by aero- 
elasticians, such as Bisplinghoff and Ashley [5.3], Fung [3.I9], 

Garrick [3 .20], Bolotin [3.4] and others. More recently, Herrman and 
Bungay [3.21] followed this usage for structural instability problems. 
When applied to nongyroscopic systems, the kinetic method is charac- 
terized by the question; What is the value of the load for which the 
general motion of the perfect system in the vicinity of equilibrium 
ceases to be bounded? This method supplies all the instabilities of 
the linear nongyroscopic system. 

When this approach is used for problems with gyroscopic forces, 
such as the rapidly rotating and deformable body, under certain condi- 
tions, one should be careful with regard to the condition that some of 
the eigenvalues are zero, A typical example is a rotating concentrated 
mass on a symmetric massless shaft, when analyzed with respect to a body 
fixed rotating coordxnate system. When X = 0 , the amplitudes are 
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bounded if the imperfections are not considered. Any sort of imperfec- 
tions with respect to body fixed rotating coordinates leads to insta- 
bility of forced vibration motion. When the same system is considered 
in a space fixed coordinate system^ there is no possibility for X to 
be zero, implying thereby, tentatively, no static divergence. However, 
experiments indicate such kinds of instability do exist [Ref. 3.9]. 

This was again attributed to the presence of imperfections (mass un- 
balance) and was explained through the familiar resonance phenomenon. 
Thus, when this kinetic method is applied to gyroscopic systems, one 
must include certain imperfections in the system. The detection of 
these instabilities associated with the imperfections is relatively 
easier in the rotating coordinate system, since they appear as X = 0 
and require additional efforts in the fixed system. It implies that it 
is always better to use a body fixed rotating coordinate system to de- 
scribe the motion. With respect to this rotating coordinate system the 
kinetic method is characterized by the question: What is the value of 

the load or angular velocity for vdiich the most general motion of the' 
imperfect system in the vicinity of steady-state equilibrium ceases to 
be bounded? Then, this method supplies all the instabilities of the 
linear gyroscopic system. At times, the solution of the characteristic 
equation may prove to be a formidable task, particularly for higher 
order systems. In such cases the Routh-Hurwitz criterion for stability 
[3.22] provides certain conditions on the coefficients of the charac- 
teristic equation. 
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The classical static approaches are the deflections of the imper- 
fection method, the nontrivial equilibrium method, and the potential 
energy method. In linear stability problems of conservative nongyro- 
scopic systems, all the critical loads are supplied by all these methods. 
With the addition of gyroscopic forces, the equilibrium method and the 
imperfection method supply the first critical load or the angular speed 
and give no indication concerning other instabilities. Actually, there 
are cases where the equilibrium may be stable although the potential 
energy is not positive definite. This effect was studied by Thomson 
and Tait [5.25], who referred to it as stability due to gyroscopic 
forces. Further, adding a little internal damping of the viscous type 
causes a loss of stability due to gyroscopic forces. Again, the sys- 
tem was shown to be stabilized under certain conditions by introducing 
external damping (air friction in the case of an unconstrained body; 
or friction at the bearings in the case of a constrained body) by 
Ziegler [5.9] and Bolotin [ 5 . 4 ], if the system is also subjected to 
circulatory (velocity independent) forces, the static methods cannot ^ 
be used. All these cases clearly reveal that static stability methods 
have shortcomings, particularly as applied to gyroscopic systems. For 
the problems considered in this dissertation, we shall assume the sys- 
tem is linear and imperfect and use the kinetic method of stability 
analysis. Whenever a simpler analysis is available, its results will 
be added to those of the kinetic method. 



5.^ DAMPING MODELS 


As mentioned in Section damping forces may be divided into 

two categories in consideration of the questions associated with sta- 
bility of the rapidly rotating systems. The first one includes forces 
caused by frictional forces from the surrounding meditnn, which can be 
considered stationary, and by contact of the rotating body with the 
adjacent fixed components. The second category includes frictional 
forces set up by the hysteresis in the material of the body, as well 
as certain other forces resulting from the relative motion between ro- 
tating parts (for example, the frictional force at the flapping hinge 
of a helicopter rotor). Damping forces belonging to the first cate- 
gory are called external damping forces and those of the second cate- 
gory internal damping forces. The internal forces must be included in 
a body fixed rotating coordinate system, and the external forces must 
be included in an inertial coordinate system. 

A precise mathematical model of damping characteristics is not 
possible due to the complex nature of the mechanism. Consequently, 
there are many models. The most common types are the linear viscous 
damping and the hysteretic or structural damping. According to the 
first model, the damping is represented by a force proportional in 
magnitude to the velocity and acting in the direction of the velocity. 
This is commonly used in vibration studies. Structural damping is 
associated with internal energy dissipation due to a hysteresis effect 
in cyclic stress . Uie energy dissipation for a cycle is assumed to be 
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proportional to the square of the amplitude and independent of the fre- 
quency. Thus, this model is valid only for harmonic motion. This 
damping produces a force proportional to the elastic restoring force 
and in phase with the velocity of oscillation, Aeroelasticians deal- 
ing with fixed wing flutter, Theoderson and Garrick [5.2?], Bisplinghoff, 
Ashley and Halfman [5.26], Fung [5.19l» Scanlan and Bosenbaum [3.27] 
among several others, commonly employed this structural damping model. 
Whirl flutter specialists who considered only external damping, Houbolt 
and Reed [3»28], Reed and Bland [5.29], used both the models; Richardson, 
McKillip and Naylor l5.30]^ Young and Lytwyn (3.31] adapted a viscous 
model. So also Bolotin [3.^] used both the models while discussing 
shaft problems with both external and internal damping included. 

Confining our attention to the rotating systems, the message from the 
corresponding references is that damping in gyroscopic systems has a 
strong influence on stability, unlike in fixed wing flutter. The com- 
parative study of Ref. 5,29 using both the models indicated that vis- 
cous damping is not as effective as structural damping in the preven- 
tion of instability. However, the overall effect of damping on the 
stability boundary does not depend on the damping model used except 
for a slight quantitative shift in the boundary. This shift was at- 
tributed to the assumption that the danqsing coefficients of these two 
models are related by equivalent energy dissipation criterion. Ac- 
cording to Ref. the relative effect of internal and external 

damping does not depend on the type of damping model used. So in the 
Investigations of this dissertation the viscous damping model is chosen 
arbitrarily. 
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Up to now we have derived the equations of motion of a rotating 
body and discussed the concept of stability and damping. In the fol- 
lowing sections we shall Illustrate the need for nonlinear elasticity 
theory, the phenomenon of frequency splitting, the advantages of using 
a rotating coordinate system, the role of internal and external damping, 
and the use of the kinetic method of stability with imperfections, by 
considering simple examples involving rotation. 

5.5 VIBRATION OF A ROTATING SHAFT 

It is well known that there seems to be a very practical limit on 
the use of a continuum model to extract meaningful conclusions partic- 
ularly in the presence of rapid rotation. However, to gain useful in- 
sight into rotational effects, some simple examples, thou^ most of 
them are not new, are considered to accomplish our stated objectives. 
The first example is a heavy uniform shaft having equal flexural rig- 
idities in the principal directions and mounted in rigid axially sym- 
metric bearings. It rotates with a constant angular velocity 
about the horizontal centerline of the bearings. Let us fix the ori-' 
gin o of the inertial coordinate system x'y'z' with the x' axis 
parallel to the center line of the bearings . The coordinate system 
xyzo rotates with the shaft with a constant angular velocity 
about ox' , This body is constrained, and hence only the elemental 
equilibrium £q. (5 .19) bo be considered. For this kind of body 
the mean axes and the principal axes coincide (see Ref. 5 •32), We 
shall fix directions oy and oz parallel to the principal or mean ' 
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axes of the shaft and initially coincident with oy' and oz" , The 
mass and elastic centers of the cross sections are assigned to lie on 
the axis of the shaft. Let v and w be the components of deflec- 
tion of the shaft parallel to the directions oy and oz , and v' 
and w' be the components parallel to oy^ and oz' (see Fig. 

The equations of motion of this shaft can be formulated either in 
an inertial coordinate system or in a body fixed coordinate system. 
First, we shall formulate them in a body fixed rotating coordinate 
system, and then transform them into an inertial coordinate system. 

For this shaft, we have the displacement vector 


7 . vf + wif - (t ^ + ng)? 

if = Q T 


( 5 . 3 ^^) 

( 3 . 35 ) 


Substituting Eqs. ( 5 * 3 ^) ( 3 . 35 ) itito Eq. ( 5 . 19 ); using beam theory 

of bending; and neglecting rotary inertia, shear deformations, and 
gravitational forces, we obtain 


m 


(tl 


dw ^ \ a V 

2 li_ ““ - vi^_ j + El — = 0 


at 


(3.36) 




m 


/ a^v sv o \ 

I — r +2 — - wfl ) 

lat^ at ° V 


2 \ 

ft - wft 1 + El — - 0 

a* 


(3.37) 
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and the following boundary conditions; 


For a pinned -pinned shaft 


v(0) = w(0) = El 


d^v(O) 

~ 2 ~ 


S^w(o) 


= El 


= 0 




vCL) = w(l) = El 


a^v(L) 


= El 




S^w(L) 


= 0 




For a pinned-free shaft 


v{0) = w(0) = El 


^^v(0) 






B^w(L) 


El 


= El 


a* 


ax 


= El 


a^w(o) 


ax 


= El 


a^v(L) 

ax^ 


= 0 


a^w(L) 


= El 




= 0 


( 5 . 58 ) 


For a cantilever shaft 


v( 0 ) - w( 0 ) = 


av(o) aw(o) 


ax 


= 0 


ax 


a^v(L) 


a^w(L) 


El 


= El 


ax 


ax 


= El 


a^v(L) a^w(L) 


ax5 


= El 


ax^ 


= 0 


We can set t) = v + Iw and write Eqs , (306)-(5.58) as 


2 4 

a n . S TJ 2 8^ 

— 5 : + (El/m) — j- - 1 ) + 2 ifi — = 0 

at ax ^t 


(5.59) 
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and the boundary conditions become: 


For a pinned-pinned shaft 


^( 0 ) 


■nCL) 








0 


0 


For a pinned-free shaft 


n(o) 


El 





For a cantilever shaft 


a^n(L) 


S Tl(0) 

El — ^ = 0 
El— ^ = 0 


Sti(0) 


a^n(L) 



Seeking a solution in the form 



(3.1tO) 


Tl(x,t) = $(x) |(t) 


bM) 
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we obtain after some algebra 


00 

v(x,t) - ^ 4^(*) cos I 
n-l * 

»„ cos [(<o^) 




+ B 


o t + 

2n 2n 


(5.^2) 


w(x,t) = V(») ( \ cos [ + I 1 

n==l ' 

where 

(“r)io - «0 * (“NR)n 

(“R>2n ' - * (“0R>n (iM) 


and ^ ®n ^ ^In ^2n constants to be determined from 

initial conditions. The coordinate functions, j and the non- 
rotating frequencies, , depend on the type of boundary con- 

ditions, and they can be found in Ref. 5.24. Equations (5.44) are 
plotted for the first mode in Fig. 5.2(a). It shows that each non- 
rotating frequency, ^®jjR^n ^ split into two frequencies, (to ) 

R In 

and (05j^)2n • This phenomenon is referred to as frequency splitting 

in Ref. 5.6. In the mode with the frequency (oi^)^^ , a point on the 

shaft traces a path in a direction opposite to that of rotation, and 
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in the other mode the point moves in oa j* 

point moves in the same direction. The former 

is called backward whirl or precessinn an^ i 

or precession and the latter forward whirl or 

precession. When > ftp \ hr>i-Vi <-v. j 

o ^ NR^n " ®odes are backward; and when 

o ^“wR^n ^'1®* (5.^2) and (3.43) show that the displacements v 
and w are bounded. Hence, the kinetic method of stability without 
imperfections fails to predict the practically observed instability at 
this condition. As stated earlier, we must introduce imperfections to 
extract all the instabilities of the gyroscopic system. These imper- 
fections, such as initial lack of straightness, will lead to some con- 
stant tenns on the right-hand sides of E,s . (3.36) and ( 3 . 37 ). Then 

% ~ the steady forced motion is unbounded. 

The same problem in inertial coordinates leads to 




m — ^ + El = 0 


^ + El — 5 — Q 


(3.45) 


(5.46) 


which are the same as those of a nonrotating shaft. The frequency is 

1 . Of . Haro .Uo, If tHa fo perfear, thara 

la »o possfMUr, fa, tpataBUfr,. Howavar, fha praaaaca of f^ar- 
factions, Uka rotating unbaianca nass, will Uad to hamonlc axciting 
torca ta^a on tha right-hand sidaa of E,a. O.h,) ihan 

tha forcad notion is unbonndad at K = (oj ) , 

hnown rasonanca phanonanon nsad to axplain tha aritical spaads . lha 



frequencies of the shaft are unaffected by the rotation. For compari 
son, the variation of frequency in any plane passing through ox' 

with is shown in Fig. 5 . 2 (b). It is clear that there is no pos- 

sibility for frequency splitting. 

Now let us introduce both internal and external damping to the 
system. We assume that these damping coefficients are equal in all 
the planes passing through ox' . If mc^ is the external viscous 
damping, the damping forces in y and z directions are given by 

/^t) and mc^(^w'/^t) . Using the transformation between 

inertia.1 and rotating coordinates 


J 




COS n t 
o 


sin n t 
o 


sin Q t cos t 

o o 


i?' 


(3M) 


the damping forces in the y and z directions are me (^ - w) 
and mc^ respectively. If me. is the internal damp- 
ing coefficient, me, ^ and tnc ^ a-ra j 

* 1 ‘^i damping forces in the 

y and z directions. Adding these additional forces to Eqs. (5.36) 
and (3.37), we obtain 




2 Sv 


2 

- 2fi — 

0 

- vfi + c. — + C 1 

° "at A 

' 0 

- U V 

at 

St 

)at ° 


+ (El/m) — jj- = 0 


(3.i+8) 
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and 


^2 

— + 2H - — - wfl + c. • — + 

o -^^ o i 



+ 


(El/m) — jj- 


0 


( 5 .^ 9 ) 


To demonstrate the role of damping we shall seek a solution in the 
form 


^ H *n(*) 

n=l 

CO 

" X] 

n=l 

where the <I>^(x) are the characteristic functions of the nonrotating • 
shaft, and and are the functions to be determined. Substi- 
tuting Eq. {5.50) into Eqs. and ( 5 .^ 9 )/ we obtain the following 

equations of Galerkin^s method: 






0 


(3.51) 
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where ~ ^ solution of Eq, (5 ■51) is given by 



where 


r + 
n 


IS 


n 


< 1 

I 


(c, + c )‘ 
1 



(3-55) 


Here r^ and are real constants, and F and G are complex con- 
stants . 

The stability of the free motion can be studied by Routh’ s method, 
and it is shown in Ref. (Page I58) that the condition for stability 
is 



(5.51^) 


The F term in Eq. (5.52) represents the backward mode, the path 
of which is an inward spiral. However, the G term represents an in- 
ward spiral motion until the condition of Eq. ( 3 . 5 !^) is met, and an 
outward one if it is not. The direction of the spiral is forward if 
^0 ^ ^"nrV backward if , Xhe same conclusions 

could be drawn if the analysis is made with respect to an inertial 
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coordinate system (see Ref, 3*9? Page 101). In both the coordinate , 
systems we need to introduce imperfections to explain static insta- 
bilities, But, as said earlier, the job is relatively easier in the 
rotating coordinate system. 

If the external damping is large compared to the internal damp- 
ing, the shaft is stable for : whereas, if the internal 

damping is very large compared to the external damping, the shaft is 
unstable for ^ , This indicates that the internal damp- 

ing has a marked influence on stability in gyroscopic systems. 

5.6 VIBRATION OF A ROTATING UNSYMMETRIC SHAFT 

In this section we shall consider an unsymraetric uniform shaft 
mounted in rigid axially symmetric bearings . The lack of rotational 
symmetry leads to time dependent coefficients in the equations of mo- 
tion if formulated in an inertial coordinate system. So, a body fixed 
rotating coordinate system is ideal for the analysis. However, the 
unsyrametry in stiffness precludes the use of complex variable tech- 
nique. This problem was studied by Ariaratnam [5 .55 3* Without going 
into the details, the principal findings that are relevant to our study 
are summarized. The internal damping plays a key role in the stability 
of the shaft. There is no question of frequency splitting since the 
shaft is not symmetric. The kinetic method of stability predicts all 
the instabilities without considering any imperfections in the system, 
unlike the case of a symmetric shaft. Finally, there is no need for 
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nonlinear theory of elasticity, since the shaft does not experience any 
steady-state equilibrium loads. 

5.T VIBRATION OF A ROTATING BEAM 

When a deformable beam rotates about an axis perpendicular to its 
own axis, it experiences equilibrium stresses in the steady-state con- 
dition. A common approach to the derivation of the vibrational equa- 
tions of motion is to treat the centrifugal loads at the steady -state 
condition as an applied load. Consequently, the need for nonlinear 
theory of elasticity was not mentioned in Refs. 5.24, 3.26, 3.54 and 
5.55 among several others. Now, we shall derive these equations using 
our general theory employing nonlinear elasticity and discuss vibration 
and stability of this beam. 

Let us assume that the beam is torsionally rigid without any ini- 
tial pretwist. We shall fix the origin o of the inertial coordinate 
system rf y' z' at one end of the beam with x' along the axis of the 
beam in the initial condition. This beam rotates with a constant angu- 
lar velocity 2^ about the z' axis. Initially, the x axis of the 
beam-fixed rotating coordinate system xyzo coincides with the x" 
axis. We shall fix the directions oy and oz parallel to the prin- 
cipal or mean axes (both are coincident for this case). It is also 
assumed that the mass center and the elastic axis lie on the axis of 



the beam. Then, we have 


r 

u 


— ^ 
q 


o 


I^x + 


T y 


+ k 2 


T> 
- 1 


/ aw av 
K — + n — 
\ ax ax 




*7^ 

+ J V 


k w 


(5.55) 


From Eq. (5.15) the equilibrium stress field is given by 


(o’ ) 

O XX 


X i 


1 

A(x) 


X 

J* m (x) dx i^ i* 
L 


(5.56) 


where A(x) is the cross-sectional area of the beam. Substitution 
of Eqs , ( 5 . 55 ) ( 5 . 5 ^) ittto Eq. ( 5 . 19) routine execution of the 

integrals after neglecting shear deformations, rotary inertia terms 
and gravitational terms produces the equations of motion for flapwise 
and chordwise vibrations as 


a 






(5.57) 


0 ( 5 . 58 ) 


and the boundary conditions become: 
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For a pinned-free beam 


v(0) = w(0) = El 


a^v(O) 


zz .2 
3x 


a^w(o) 


= El 




= 0 


a^v(L) 


S^w(L) 


El 


zz .2 


= El 


yy . 2 




= El 


zz 




( 3 . 59 ) 


S^w(L) 


= El 




= 0 


For a cantilever beam 


v(0) = w(0) = 


^v(O) aw(0) 


= 0 






S^v(L) 


El 


zz 2 


= El 


a^w(L) 




= El 


a^v(L) 




= El 


yy 


a 

“ 


- 0 


(5.6o) 


where 



( 5 . 61 ) 


El and El are the flexural rigidities of the beam in two prin~ 
yy zz _ ^ ' 

cipal planea; m is the mass per unit length; and L is the length 
of the beam. The Eqs . (3 «5T) - (3 .^l ) are self-adjoint and they are un- 
coupled. Consequently, we can find normal coordinates in each principal 
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plane. No coriolis forces are present, and hence the problem is non- 
gyroscopic. The centrifugal inertia forces enter only in the equation 
representing vibration in the oxy plane. As a result, the frequen- 
cies with rotation in two principal planes, even for a symmetric circu- 
lar beam, differ, being higher in the flapping plane ozx and lower in 
the edgewise plane oxy . This can be interpreted as frequency split- 
ting due to centrifugal inertia forces, if the beam is symmetric. How- 
ever, there is no possibility for instability of this elastic beam. 

But, when rigid body motions [v = xp(t) , w = xq(t)] for a pinned- 
free beam are considered, from Eqs. ^nd (5.58) it is seen that 

the frequencies with rotation are (60„) = and fe_) =0 , which 

clearly indicates that the beam is statically unstable in the chordwise 
direction. Loosely, this instability is analogous to that of the sym- 
metric shaft when = (a>„) . It can be seen that the kinetic method 

of stability without imperfections can predict the present instability 
because the system in nongyroscopic. Though the system of Eqs, (5 •57)- 
(3.61) is self-adjoint^ no exact solutions are possible, even for a 
uniform syimetric beam. However, approximate solutions by Galerkin 
and Raleigh-Ritz methods are given in Ref. 3 0^ and Numerically, 

exact solutions will be given in Chapter IV for a cantilever beam. 

3.8 VIBRATION AND STABILITY OF A ROTATING BEAM WITH CONING ANGLE 

The problems discussed in Sections 3*5 "5*6 turn out to be spe- 
cial cases of a broader problem, whose analysis throws more light on 
both of them and elucidates the nature of the coupling in such systems 
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generally. This model, as shown in Fig, consists of a slender 

beam capable of bending in two planes perpendicular to its neutral 
axis and spinning about an inertial fixed z* axis, which makes an 
arbitrary angle (^ " Pq) with the axis of the beam, is called 

the coning angle of the beam. Due to this angle, there enter coriolis, 
centrifugal and steady-state equilibrium forces into the equations of 
vibration. From the academic point of view, this problem gives an op- 
portunity to study the effect of all these forces on vibration and 
stability; and from the practical point of view, these results can be 
profitably used for a better understanding of whirl flutter in prop- 
rotors, Tdierein the blades of certain rotors have built-in precone in 
order to reduce the out-of-plane bending moments and feathering bear- 
ing loads . 

Referring to Fig. 3*3^ we have 

7 = r (- 11 — - 5 — \ + (3.62) 

V ^x / 

^ il sin 6 + cos P ( 3 »^ 3 ) ■ 

O o *^0 O O v-' / 

L 

cr^ = [1 /a(x)1 J* cos^ ni(x) dx ? i^ (3*^^) 

o 

where A(x) is the area of cross section of the beam. We have neg- 
lected longitudinal deformation in addition to shear deformations . 
Substituting Eqs . (5 .62) -(5 . 63 ) into Eq, (3.19)# after neglecting 
rotary inertia and assuming that the rotations are larger than the 
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FIG* 5 •5*"‘'Rc)tating beam Inclined at an arbitrary 
angle with the axis of rotation. 


- 5 ^ - 


strains^ we obtain the following: 


m 


/3^v 

^2\ 

r.2 2 

av 

( 2 ’ 

2 — sin p - vH I 

- cos p 

0 0 

T(x) — 

vat 

at / 

ax 


^ ^ ^ \ 

-2 "^zz-2 = " 

Sx \ / 


(5.65) 


m 


/ Sv 

I — 2 + 2 — 

\at^ at 


2 2 

sin 6 - wH sin p 

000 


) 


22° 

- n“^ cos p — 
o o . 

ax 


av 

T(x) — 

ax 


•?KS) 


X fi m sin ^ cos p 
o o o 


(5.66) 


where 


J-i 

T(x) = / m(x) xdx 


(3.67) 


The boundary conditions for a plnned-free and a cantilever beam are 
given by Eqs . (3*59) and (5*60), respectively. 

The above Eqs. (3*65) and (3*66) and the appropriate boundary 
conditions can also be obtained by equilibrium considerations for an 
element of the beam considering the equilibrium centrifugal load as 
an applied one. The present method explicitly demonstrates the need 
for nonlinear theory of elasticity when the system is subjected to 
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equilibrium steady-state loads. Now^ one can safely conclude that any 
general theory of vibration of a rotating deformable body must include 
nonlinear theory of elasticity. In some simple cases, for example, in 
the case of a rotating shaft, the terms associated with the nonlinear 
theory may be dropped, since it experiences no steady-state centrifugal 
loads. The right-hand side in Eq. (5-66) is constant, and hence for 
the study of vibration this can be dropped. An exact solution of these 
equations in closed form appears to be not possible. However, we shall 
concentrate on approximate solutions in this section and shall present 
a numerical solution in Chapter IV. 

First, we shall show that the approximate frequencies of a uniform 
coned beam with rotation can be deduced from the information on approxi- 
mate frequencies of the same beam without coning presented in Ref. 
Second, the stability of the uniform symmetric coned beam for a range 

of values of B and Q will be discussed. We assume a solution of 
o o 

Eqs. (3.65), (3.66), (3.59), and (3.60) of the form 




( 5 . 68 ) 


n=l 




(3.69) 


n=l 


where characteristic functions of the self-adjoint 

boundary-value problem 


El ^ 




zz 
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and the boundary conditions become; 


For a pinned-free beam 


^^(0) 


a\(L) 


s\(0) 

(L) 


ax 


ax5 


= 0 


=: 0 


For a cantilever beam 


( 3 . 71 ) 


$^(0) 


a\(L) 


ax 


a\(o) 

ax 

a^\(L) 

■y" 


= 0 


= 0 


We then obtain 


i’n - 


RT]^n o o o 


p =0 


q ‘h 26 0 sin p + 
^n o o 


sin^ p 

' R^^n o o 


I-- 


0 


(5.72) 


where the the frequencies of the corresponding boundary- 

value problems (5*70) and (5.71) with EI^ distribution instead of 
Elg 2 • Though the frequencies are different, the characteristic func- 
tions are the same. The values (a^jj^)jj and )^^ , which are equal 
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if the beam is symmetric^, can be obtained from the charts of Ref* 5 ,5^ 


just by considering the rotation as cos instead of * It 

can be seen from Eqs. that, even for a symmetric beam, we cannot 

use the complex variable technique to solve the equations as in Section 
5*5* Seeking a solution of the form 


p P e 


Xt 


= Q e 


Xt 


(5.T5) 


where X ^ p ^ Q are constants, we find that the exponent X is given 
by the characteristic equation 


X + 


sin^ 0 

' RT] 'n •^n o o 


(^T, 


“o (“EC'n ■ "o Po 


= 0 


(5.TU) 


which leads to 




1 

2 




± { sin^ g cos^ p + (03„.)^ 

' o ^ RT)-'n o ' R^'^n 




- KOl " 


2)i 


(3.75) 
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p p 

The general solution is bounded if (>-i ) , ) are both 

^ In'^ ’ '■ 2n'^ 

negative. This condition is satisfied when sin B lies outside 

0 o 

the intervals 


K f 

^ RT]^n 


I cos B 
o 


2 2 

sin P 
o 




RC^n 


(3.76) 


This inequality imposes certain restrictions on ^ and 0 for 

o o 

stability to prevail. Let (p ) be the minimum value of B for 

' o'vn 

p p 2 2 

which sin , and the corresponding value 

of n be } snd let and (Q ) be the corresponding 

o o vn o wn o wn 

2 2 2 

values at which (^g^^ illustrate the stability we 

consider a viniform symmetric beam and also divide the interval 0 to 

7 t /2 of into several subintervals . Further we assume the beam is 

such that 0^)^^ < 0^)^ < < ... , and < <‘^o)y 2 < ••• • 

with the same inequalities assumed valid for (6 ) and (J2 ) . 'The 

^ o'vm ' o'^wn 


first subinterval of 6 is 0 s B < fB I , 

o ^o '•^o'^vl 


In this interval both 


'•*'"^RT|^n ■ ^o ^o^ greater than sin^ , and 

hence the beam is stable for all the values of H . At B = (B ) 

o ^o ^^o'^vl ^ 

and = ^^o^vl eigenvalue X = 0 is a double root of Eq. (3,7h). 

This corresponds to a solution of the form 


Pi = Pi + Pgt 

'll = \ (5 .77) 
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which shows p. is unbounded. With (P ) , constant and > (J2 ) . , 

^ o'^vl o ' oVl ' 

the value of ft sin 6 lies inside the interval given by Eq, (5.76), 
o o 

and hence the beam is unstable. This is illustrated in Fig, 5*^(a). 

In the second subinterval (P ) , < P <(?>)■, , an increase in 

^*^o^vl o ' o'^wl ’ 

from reduces the angular velocity ft^ at which instability 

sets in first and the bandwidth of instability continues to be infinity 

until P = (6 ) _ and ft = (ft ) , at which point the solution is 
o ^ o'wl o ^ o'wl ^ 




= Qn 


+ Q2t 


(3.78) 


It shows q. is unbounded. With further increase in ft and 
1 o 

P = (P ) 1 > the value of ft sin p lies outside the interval, 

o ' o'wl o o 

Eq. ( 5 . 76 ), and hence the beam is stable. Thus, for P^ = (P^)^^ 

the infinite bandwidth of instability of the mode begins to change 

to a finite one. The variations of frequencies of the first two modes 

with ft^ for this range of p^ are shown in Fig. 5«^(h). 

Next, subinterval (p ) , < p < (p ) . is considered. For each 

' o'wl o ' o'v 2 

P , there are two values of ft : namely, (ft ),, and (ft )^, 

o' o ' o'll ' o '21 

between which the shaft is unstable. The range between (ft ),, and 

o' 11 

(ft^)gl is called the bandwidth of instability of the first mode. An 

increase in P^ reduces the bandwidth of instability of the first mode 

and the beam is stable for all speeds ft > (ft . When P = (P ) 

o'' o '^21 *^o '^o‘^v 2 

and ^^o^v 2 ^ beam is again unstable statically at 

fto = (^q)v 2 ’ beam has finite bandwidth of instability in 
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the first node and an Infinite one in the second mode, as shown in Fig, 

In the fourth subinterval (P ) < p < (p ) _ , the first mode 

o^v2 o o w2 

has a finite bandwidth of instability^ whereas the second mode has an 
infinite bandwidth until p = which the infinite bandwidth 

of the second mode begins to change to a finite one. At this point for 

any beam is stable. This process repeats as we in- 

crease p . Eventually, at a given P and for all values of n >n , 
o o o' 

where n is the value of n for which (fl ). , S (Q , these 

o ' o^ln +1 o'2n ' 

o o 

bandwidths of instability overlap ones immediately preceding them, and 

hence for all the speeds S: the beam is unstable. This is 

*0 

shown in Fig. 

Finally, when p^ = tt/ 2 , these bandwidths shrink to points when- 
ever which is the familiar shaft prob- 

Im. A. pointed out earlier, the anvlltudea p^ and are bounded 

if imperfections are ignored. This case is shown in Fig. 5.4(f). 

The next problem is to define the symmetric beam model by consider- 
ing rotary inertia and shear deformations . This problem with = 0 
was studied by Pederson (Ref. 5.57), who also introduced both internal 
and external damping, and anisotropic elastic restraints at the bear- 
ings. Here, there is no need for nonlinear theory of elasticity, since 
the equilibrium stresses are absent. A space fixed coordinate system 
was used because of the anisotropic bearings. Therefore, the insta- 
bilities of the system were brought out through the resonance phenom- 
enon by considering imperfections. Instead of going through all the 
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details, we shall sunmiarize the principal findings that are relevant to 
our study. 

When rotation speed is zero, the system frequencies can be di- 

vided into two groups. One corresponds to the slender beam group, and 
the other to the rotor group . These two groups are coupled even in the 
absence of rotation. Therefore, the system has an infinite number of 
frequencies. When this symmetric shaft with isotropic bearing restraint 
is spun at a speed , each frequency is split into two frequencies; 

the higher one with respect to inertial observer is called forward vAiirl; 
and the lower one, backward whirl. Thus, at any speed , there are 

an infinite number of forward and an infinite number of backward fre- 
quencies. So, the frequency splitting takes place even when the system 
is studied with respect to an inertial coordinate system. This split- 
ting is due to gyroscopic moments. If the imperfections are ignored, 
the system is stable at all speeds vdiich is a typical characteristic of 
gyroscopic systems under certain conditions. 

When forced vibration of the shaft due to static unbalance was con- 
sidered through modal analysis, it was shown that if the bearings of the 
rotor have anisotropy, the amplitudes of both forward and backward whirl 
modes are unbounded. Aiwther interesting observation is the role of 
both internal and external damping. With the introduction of these 
dampings (intemal/extemal = 0.2), the amplitudes of the backward 
whirl modes are reduced more than those corresponding to forward idiirl 
modes. However, for the isometric bearing case the amplitudes of the 
backward whirl modes are always bounded both with and without damping. 
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5.9 OTHER RELATED EXAMPLES 


The other problems of Interest with a continuum mathematical model 
are the transverse vibrations of a spinning circular disk. Refs. 5-58 
and 5*39> plane vibrations of a thin uniform spinning ring. Refs, 
5 .^ and 5.^1; and the stability of a spinning body containing elastic 
parts. Refs. 5-18^ 3.^1^ and 5*^2. In the disk problem, the motions of 
all its points are perpendicular to the plane of the disk so far as 
first-order terms are concerned, equilibrium centrifugal loads are pre- 
sent, and hence there will be some increase in vibration frequencies 
with rotational speed. But there is no frequency splitting as there 
are neither centrifugal inertia nor coriolis forces . The disk is 
stable for all speeds. In the ring problem, the equilibrium centri- 
fugal loads increase the frequencies, and the coriolis forces cause the 
nodes not to remain fixed in the ring. The effects of rotation on the 
third problem depend on the orientation of elastic members in the body, 
A comparative study of stability methods based on the Liapunov direct 
method is described in Ref. 5*^3. 

3.10 summary 

For a general theory of vibration of a rapidly rotating deformable 
body, it is always safe to include the nonlinear theory of elasticity. 
Though this fact was not explicitly mentioned in the literature dealing 
with rotation, the additional terms in the equations due to nonlinear 
theory are brought in through an effective applied load concept which 
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cannot be easily Implemented for a general elastic body. It is clear 
from the examples that the frequency splitting phenomenon depends on 
the axis system considered, the geometry, and the orientation of the 
body with respect to the rotation vector. In general for this phenom- 
enon to take place, the body must have different modes with equal fre- 
quencies in the absence of rotation, and the rotation must add either 
centrifugal or coriolis terms to the equations of motion. The boundary- 
value problem governing the vibrational motion about the steady-state 
condition is not in general self-adjoint, and hence the modes are not 
orthogonal. With regard to stability criteria, the kinetic method ^en 
applied to a rotating elastic body falls in some cases if imperfections 
are not considered. Therefore, it is always safe to include imperfec- 
tions in the system and apply the kinetic method. Another stability 
method, if possible, is the Liapttnov direct method. The instabilities 
associated with Imperfections are relatively easier to detect in the 
body fixed rotating coordinate system, if one is suitable, rather than 
in the inertial coordinate system since they appear as X = 0 in the 
rotating system. The last feature is the role of internal and external 
damping. It was seen that internal damping, when it exceeds external 
damping, has a destabilizing effect on stability in gyroscopic systems 
under certain conditions. Finally, the existence of several interesting 
bands of instabilities of a rotating beam with coning angle as a param- 
eter is established. This analysis is believed to be new. Chapter VI 
will examine the effects of coning of the blades of a prop-rotor, and 
of internal damping at the blade flapping hinges on whirl flutter 
stability. 
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IV. NUMERICAL CALCULATION OP BENDING FREQUENCIES 
OF A ROTATING BEAM WITH CONING ANGLE 


4.1 INTRODUCTION 

The governing equations of bending vibration of a rotating beam 
with coning angle were derived in Section 5.7^ and it was also shown 
that the approximate frequencies can be deduced from the information 
presented in Refs. 4.1 and 4.2 for a class of rotating beams without 
coning. In this chapter we shall present an exact method (within the 
accuracy of the numerical calculations) to calculate the frequencies 
of the coned beam with rotation, and these values will be compared 
with those of the approximate Galerkin*s method. 

4.2 THEORETICAL CONSIDERATIONS 

It can be seen fromEqs. ( 5 . 65 ), (5.66), ( 5 . 59 ), and ( 5 . 60 ) that 
the eigenvalue problem governing the vibration of the beam considered 
is a two-point boundary value problem since all the boundary values 
at the root are not known. So a solution of these equations for a 
set of arbitrarily chosen unknown boundary conditions and the fre- 
quency are first obtained at a given rotational speed by the Runge- 
Kutta numerical integration process [4.5) after transformation of 
these equations into a set of simultaneous linear first-order 
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differential equations. Ihe actual solution of the equations is a 
linear combination of all these arbitrary solutions. So imposing the 
condition that the actual solution must satisfy the known boundary 
conditions at the tip of a cantilever beam results in a set of homo- 
geneous equations. The vanishing of the determinant of their coef- 
ficients gives the frequency at the preselected rotational speed, and 
the mode shapes along the beam can be obtained from the linear combi- 
nation of these solutions at this frequency. The steps involved in 
this process are as follows: 

(a) Transformation of Eqs . ( 5 . 65 ) and (3.66) into state 
vector form; these are given by 
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and 


>9 

w(x,t) 

v(x,t) 




M 



yj^(x) sin OJ^t 
y^(x) cos 

V(%r)j 

U I , h 
c El / (m L ) 

yyo o ' 

I /l 

yyo' zzo 


(^.2) 


The parameters m , , etc , are no ndlmens localized with respect to 

their corresponding values at the origin, and the variable x is non- 
dimens ionali zed with the length L of the beam. The boundary condi- 
tions are given by 


^1 = ^2 = ^5 = y6 = ° X = 0 (b.5) 


Vj, Yk = ^ Vq = = 0 at X = 1 (b.lf) 


- 72 - 



(b) Values of the rotational parameter and the frequency 

parameter v„ are assumed. 

(c) Starting at x = 1 , the uncoupled differential equation 

y' = - mx is integrated by the Runge-Kutta step-by-step 
numerical Integration process. 

(d) The four known boundary conditions at the root are set 
to their respective values, and the remaining boundary 
conditions are assigned arbitrarily as 


set 1 

^3 = 

1 , 

yi. - y? * ys ■ ° 


set 2 

^4 = 

1 , 

y3 ■ y? ' yg = ° 


set 5 

II 

1 , 

y3 ■ yii ' ya ' ° 


set 4 

^8 = 

1 , 

yj - yi, ' yy = 0 

(4.5) 


Starting from x = 0 , four sets of solutions of Eqs. (4.1) are ob- 

tained by the Runge-Kutta numerical integration process. At the tip, 
we shall write 

y. .(1) = E(i/j) 

^7 J 

i = 1,2,... ,8 

j = 1,2,5, 4 (4.6) 
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where suffix j indicates the solution of the equations with the 
jth set of boundary conditions [Eq. (h. 5 )]. 

Since Eqs, (^.1) are linear, the actual solutions at the tip 
are given by 

h 

yi(l) = ^ (1^.7) 

j=l 

So also 

k 

>' 1,1 

j-i 

where the a^ (j = 1 , 2 , 5 , h) are constant coefficients. 

The known boundary conditions at x = 1 represented by Eq. (U.h) 
require that 

k 

E(i,j) = 0 

j=i 

(i = 5,6,T,8) (k. 9 ) 


and hence we must have 


|E(i,j)J = 0 
(i = 5, 6 , 7, 8 ) 

(j = 1 , 2, 5 , 4) (k.io) 
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The value of the determinant is zero if and only if the assumed 
frequency V- is a correct one. So the value of the determinant is 
calculated for a range of values of at a preselected value of 

V, , and it is plotted against v„ • The points of intersection on 
the frequency axis give the frequencies of the system. For each value 
of ^ there are two values of ^ 

mode because of the frequency splitting due to rotation. Both of these 
frequencies can be extracted from Eqs. (If. 2). However^ to determine 
the type of whirl, we must also solve Eqs, (^.1) by writing 


w = yj(x) cos 03j^t (^•11) 

V = y^(x) sin ojj^t (^*12) 


Then we obtain the same frequencies with different mode shapes. Guided 
by the experience from Eqs. (5.J+2) and (jAj), if we write the solu- 
tions of Eqs, (5»65)“(3«68) in the form 




n=l 


+ c cos t + d sin (n>_). t 

n ' R''2n n ^ R’'2n 


(^. 13 ) 


w = 


CO 

f ^n 

n=l 


cos (cD ) t + H sin (tn„), t 
' R‘'ln n R-'ln 


+ C cos + D sin (n>_)^t 

n ' R'^2n n ^ R^^2n 


(it.iit) 
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the coefficients ^ determine the direction of 

whirl of the mode with frequency , and the coefficients c , 

R Xn n. 

d , C and D determine the direction of the other mode. If the 
n ' n n 

beam is symmetric identical. 

Knowing the frequency, the coefficients > a^ and ai^ are 
determined from Eqs. (^.9) for an arbitrary value of a^^ , say 1 , 

and the mode shapes are obtained from Eq. (4.8). 

These equations are programmed on an IBM 560 / 67 , and the exact 
frequencies and mode shapes of a uniform, symmetric beam are obtained 
for different values of the parameters and 

The assumption of the parameter v,, in step (b) above is guided 
by an approximate one-term Galerkin solution. In this method, the fre 
quencies of each mode of a uniform beam are calculated by using the 
corresponding nonrotating mode shape. These frequencies are also in- 
cluded in Figs, 4,1 and 4.2, 

4.3 RESULTS AND DISCUSSION 

Both the exact and the approximate one -term Galerkin solutions 
are plotted in Figs. 4. 1-4 ,3 for a uniform symmetric beam. Also in- 
cluded in Fig. 4.1(a) is a five-term Galerkin solution from Ref. 4,1. 
It is seen that the approximate one -term results are in good agreement 
with the exact results for the second mode. For the first mode, how- 
ever, the maximum error with one mode is larger, about in the fre- 
quency, and that with the five mode is less than 3 ?^. From these re- 
sults it appears that a one-term Galerkin method can be used for a 
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. k,2 — Influence of rotation on the first bending frequency of a 
uniform beam with as parameter. 


FIG 





rapid estimation of the frequencies of a rotating coned beam. 

At any , the effect of rotation is to split each frequency of 

a nonrotating symmetric beam into two and alter them with • In 

fact, the rotational parameter v-y , depending on the value of , 

may cause one of the frequency parameters zero or 

imaginary. If (^p^)2n imaginary, then the beam is dynamically un- 
stable; and if it is zero, the beam is statically unstable. Both of 
these instabilities associated with the nonrotating first mode are 
shown in Fig. with as a parameter. It is clear from the 

figure that the type of instability or stability, as discussed in Sec- 
tion 3.7, depends on the range of . 

The effect of rotation at different values of on the vibration 

mode shapes is shown in Fig. 4.5. An examination of these shapes indi- 
cates that, although some differences between the modes exist, they are 
relatively small, particularly for the second mode. We shall use the 
exact flapping frequencies of the rotating uniform beam in the follow- 
ing chapter. 
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V. LEAST WEIGHT DESIGN OF A ROTATING CANTILEVER 
BEAM WITH ONE FIXED FLAPPING FREQUENCY 

5.1 INTRODUCTION 

The main object of structural optimization is to achieve minimum 
weight designs subjected to certain design requirements. These re- 
quirements may be termed constraints on the optimization process. More 
emphasis is placed on the minimization of weight because of the over- 
vdielming importance of this parameter, particularly in aerospace ap- 
plications . 

For optimization studies there are two mathematical models for 
representing the structure. One is the elastic continuum model, and 
the other is the distributed mass finite element model. The former one 
is suitable for simple, continuous structures and the latter for all 
types of structures. Once the problem, the type of constraint, and the 
mathematical model of the structure are known, the next step is to find 
a suitable optimization technique which accomplishes our stated objec- 
tive, In the literature there exist several kinds of such techniques 
depending on the problem and the tastes of the investigator. 

When the mathematical model of the structure is of a continuum 
type with differential equations as constraint, one method of optimi- 
zation is to use control theory techniques, which are themselves 
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ramifications of the calculus of variations. The constraint equation 
is an ordinary differential equation with one independent spatial var- 
iable if the structure is one -dimensional like a beam^ and it is a 
partial differential equation for a two-dimensional structure such as 
a plate. Many simple problems can be described by differential equa- 
tions, particularly in the field of optimal control theory. In the 
case of a discrete model, mathematical programming techniques have 
been used. Both of these approaches have been developed, and they 
represent two complementary methods of solution. 

model is adapted for the problem of interest in this 
chapter. This model was widely used by earlier investigators, Turner 
[5.1], Keller and Niordson [5.2], Niordson [ 5 . 53 , Olhoff [ 5 . 4 ], Sheu 
and Prager [ 5 » 5 ]» Prager and Taylor [5.6], Ashley and Mclntosy [5.T], 
Armand and Vitte [5.8], Weisshaar I5.9], Armand [5.10] and Vepa [5. 11], 
to name only a few. The literature on the second model, which is not 
our concern here, can be found in Refs. 5 . 12 to 5.I6 among others. 

In this chapter the constraint equations and boundary conditions 
of the stated problem are deduced from those of the coning beam of 
Section 3.7, and then a set of nonlinear differential equations and 
boundary conditions, which are necessary for the optimal depth distribu- 
tion, are derived using optimal control theory techniques. Finally, a 
numerical solution of these equations is obtained using the ’’transi- 
tion matrix" method described by Bryson and Ho [ 5 . 17], and a physical 
interpretation of the optimality condition is given. 
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5.2 STATEMENT OF THE PROBLEM 

The problem considered is a beam rotating at a constant speed with 
its axis perpendicular to the axis of rotation. The beam has a fixed 
width and length with a solid rectangular cross-section. Our objective 
is to find a distribution of depth along the span such that the weight 
is optimum and the fundamental flapping frequency with rotation is the 
same as that of the uniform depth beam with the same rotational speed 
and geometry. Two sets of boundary conditions can be considered for 
this rotating beam. One set of boundary conditions will be those due 
to clamping at one end of the beam while allowing the other end to be 
free. The second set of boundary conditions will be due to fixing the 
beam with freedom to rotate at one end while allowing the other end to 
be free. This set of end conditions is referred to as "pinned-f ree" , 
However, the optimum depth distributions will be obtained only for 
clamped-free boundary conditions. 

The equilibrium equations governing the flapping vibrational mo- 
tion are given by Eqs . (3.57), (3.59), (3.60), and (3.6I) with the 
chordwise elastic displacement v being zero. Setting 

w(x,t) = w(x) e ® 

Eqs. (5.57), (5.59), (5.60), and (3.61) can be written in nondlmensional 
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form as 


[t^(x) w‘^^(x)]'^ + [xt(x) w' (x) - w^'' (x) u(x')] 

P li 

- C t(x) w(x) = 0 

u'(x) = - xt(x) 

where 

w(x) = w/l 
X = x/l 

t(x) = T^(x)/t^ = local depth/depth at the root 
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A closed form analytical solution of the boundary-value problem 
[Eqs. (5*2), (5*3)j and (5.5)3 l-S not possible even for a uniform beam. 
However, the frequency ^ selected is obtained from 

Fig. 1^.1. 

Now we are in a position to state the optimization problem. We 
want to minimize the integral 

1 

J = y* t(x)dx ( 5 . 6 ) 

0 

representing the dimensionless mass of the beam. J is also called a 
performance index or merit function. The depth function is subjected 
to the additional conditions represented by Eqs. (5.2), (5.5), and ( 5 . 5 ). 
This is a problem of the calculus of variations and may be viewed as one 
of the Bolza type (Ref. 5 .I 8 ). This problem can be reduced to a more 
conventional type by converting Eq. (5.2) into first-order equations. 

Let us introduce the new variables 

p = w' 

q = t^w**^ 

r = (tV')" ( 5 . 7 ) 
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Then 



p' = q/t^ 

u' = - xt (5.8) 

with boundary conditions 

w(0) = p(0) = r(l) = u(l) = q(l) = 0 (5.9) 

The problem represented by Eqs. (5.6), ( 5 . 8 ), and (5.9) is a kind 
often encountered in optimal control theory (Refs, 5 . 17 ). However, 
here the role of time is being played by the dimensionless coordinate 
X . The control function here is t(x) , and the state variables 
are r , q , w , p , and u . 

5.3 REVIEW OF OPTIMAL CONTROL THEORY TECHNIQUES AS APPLIED TO THE 
PROBLEM CONSIDERED 

The general problem is to choose a control function t(x) and 
the state variables yj^(x) , yg(x) , ... , y^^(x) which satisfy the 
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conditions 


?(X) 

= f(t,y,x) 



yjo) 

= 0 



y,(i) 

= 0 

J “ m "t" \ y 

(5.10) 


and which also minimize 


1 

3 = J" t(x)dx (5-11) 

0 

1 

The function f is required to be differentiable with respect to y 
and t y but may be piecewise continuous with respect to x . The 
functions y are required to be piecewise continuously differentiable^ 
but t may be only piecewise continuous • 

The classical approach to the solution of the problem is to ob- 
tain a list of additional conditions which the solution of the problem 
must satisfy and then construct a set of functions which satisfy all 
these necessary conditions. Since all the functions satisfy only 
necessary conditions, they are merely candidates for an extremal, 
rather than extremals themselves. So the one among the functions 
lAich yields the smallest value of J is the solution of the problem. 
The lack of sufficient conditions for this kind of problem forces the 
designer to make decisions on the basis of his experience and engi- 
neering knowledge rather than rigorous mathematics. However, the 
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distinction between necessary and sufficient conditions is not of great 
importance in problems of classical mechanics where equilibrium problems 
require stationarity of the functional rather than its extremal « Some 

’ I 

discussion on. the necessity of certain conditions for optimization prob- 
lems were discussed by Caratheodory [ 5 . 19]. 

A very general principle, which furnishes the tools for the above 
optimal design problem, was hypothesized by Pontryagin [ 5 . 20] on the 
basis of the work by him and his associates. This principle is called, 
the "maximum principle" . In this principle Pontryagin introduced the 
Hamiltonian in a more direct manner as opposed to an indirect manner 
in classical mechanics as a Legendre transformation of the Lagrangian 
density. The Hamiltonian is defined as 

H[y(x),t(x),x] = t(x) + >^{x)fEy(x),t(x),x] (5.12) 

where X is a vector of Lagrange multipliers. According to this prin- 
ciple, to find the control function t(x) that gives a stationary value 
of J , the following set of additional equations constitute the nec- 
essary conditions for an extremal: 

X' = 

St 


hy) 


0 


( 5 . 15 ) 


(5.1^) 


- 91 - 



and 


x^(l) =5 0 i = 

X^(0 ) = 0 j=m+l,...,n (5.15) 

The multipliers X are continuous and have piecewise continuous de- 
rivatives . These are also called adjoint variables, and the condi- 
tions, Eqs . (5#15)> also often referred to as transversality con- 
ditions. The Hamiltonian function H serves only as a device for 
conveniently displaying necessary conditions for the Bolza problem. 

Out of all the possible controls, one which yields the maximum H is 
the optimum one. 

5.^ APPLICATION OF OPTIMAL CONTROL THEORY TO THE PROBLEM CONSIDERED 
For the problem of Section 5.2, the Hamiltonian function 

H = t(x) - X^t(x)x + + X^q/t^ + X^r 

+ X^[- c\^(ptx - uq/t^) + (Vjj)jj c\w] (5-I6) 

leads to the following necessary conditions for an extremal; 

2 h- 

^r ^ (Cont'd) 
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X' = - 
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and 


u 


X' = 

q 


X' = 

V 


X /t^ - u/t^ 

p' 1 r ' 


(5.17) 


t = 


5>-pq + 3^yc\^qu 

1 - X^x - X^c [v^px - (Vj,)^]^ w] 


lA 


(5.18) 


\(0) = X (1) = X^(l) = X (0) = X^(o) = 0 (5.19) 


Finally, we end up with a system of ten nonlinear, first-order differ- 
ential equations represented by Eqs . (5.8) and (5.17), optimality con- 
dition, Eq. (5.18), and the boundary conditions, Eqs. ( 5 . 9 ) (5*19), 

for eleven unknowns t , q , w , p , u , X^ , X^ , X^ , X^ , X^ , and 
t . Since half the boundary conditions at the root x = 0 and the 
other half at the tip x = 1 are specified, this problem is referred 
to as a two-point boundary value problem. An examination of these 
equations reveals no possibility for an analytical solution and hence 
a numerical solution is given in the next section. 
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5.5 NUMERICAL SOLUTION 


The transition matrix approach developed and described by Bryson 
[ 5 . 17 ] has been successfully applied for some problems of this nature 
by the earlier investigators Armand [5.8] and Weis shaar [5.9]. Since 
this method has become a common approach^ we will not review it here. 
However, we shall present the steps involved in applying this method 
to our particular problem. 

(a) To obtain the transition matrix, the perturbed equa- 
tions of Eqs . ( 5 . 8 ) ( 5 * 17 ) are required, and these 

are given by 


(&u)' 

= - x&t 

(6X )' 
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p il 

- - + xvfc (t6X X 5t) 

w 1 ' r r ^ 

(&x )■' 

= - (t&x^ + X^5t) 

(8q)' 

~ Sr 

(8r)^ 

= (8uq/t^ + 6qu/t^ - 55tuq/t^ 


P ll 

- 6pxt - xp6t) + c (tSw + w6t) 

(6w)' 

= gp 


= - (8X^q/t^ + X^5q/t^ - 5Xj.q5t/t^) 


(Cont' d) 
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and 


(^q)' = - (5X^u/t^ 

+ X^6u/t^ - 3^j.uS't/t^) 

(gX )' = - 5X 

^ r' q 

(d^)' rz - 5 qStA**^ + ( 5 . 20 ) 

1 2 !)• 

[&x^q + 

+ X^q&u)] [1 - X^x - "k^c^ [v^px 

■ Apq + >^^V^c\q) [xSX^ 

+ 5X^C^ [xVjP - (Vj^)^jW] + X^xv^c^Sp 

■ \ (''r)ii 

4- |(*^t:^/3) {1 - tvjpx 

- (v^)ii»H^| 

(b) We express all the Eqs . ( 5 . 8 ), (5.17), and (5.20) in a 
compact and convenient state vector form as 


y' = £ (y, X, t ) 


( 5 . 21 ) 
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and 


vdiere 


yj^ = {Oj at x = l 

yg = {Oj at X = 0 

y^ = fO] at X = 1 

= {O] at X = 0 


( 5 . 22 ) 





(uX X qr)' 


(wpX XX)’ 
^ u q 



(5.23) 


(c) We now guess the unspecified conditions yj^(O) 

(d) Integrate Eqs. (5.21) for 5 different values of y^(0) , 
each time taking one of the perturbation quantities of 
y^(0) as unity and the other values of y^(0) as zero, 
with the values y^(o) , yg(0) , and yj^(O) being un- 
changed, Each time one column of the transition matrix 

T (5 X 5) is filled up by the values of y^ at x = 1 . 
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(e) Then, calculate a change Ay^ in at x = 0 such 

rw 

that the error in y^^ at x = 1 is reduced by an amount 
e ; it is given by 


r'J 

= - € T 

(at X = O) (at X = 1) (5.24) 


(f) Calculate the new initial conditions 


?! = Vi + (5.25) 

(new) (old) 

other values being the same. The same process is re- 
peated until y^ = [o] at x = 1 to the required ac- 
curacy . 

The main problem in this iterative integration process is to guess 
the initial values for a rotating beam. Presently, this is overcome by 
a sort of perturbation on the rotational frequency parameter. First, 
the optimum solution of a nonrotating cantilever beam is obtained by 
estimating the required initial values from the eigenfunction solutions 
of a uniform depth beam. Since the eigenvalue problem of the nonrota- 
ting cantilever beam is self-adjoint, the adjoint variables are linear 
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functions of the state variables . That is. 
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(5.27) 

and the boundary conditions associated with these variables lead to 


(5.26) 



(5.28) 
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where A is an arbitrary constant. The above result was reasoned by 
Tuimer [5.I]. Making use of this relation and the eigenfunctions of 
the nonrotating beam, the ratio of r(o)/q(0) = -1.5775^ > was esti- 

mated by Weisshaar [5.9]. An estimate of q( 0 ) for a uniform depth 
beam can be obtained from Eq. (5. 18 ) as 

t^(0) - 3 \(o) q(o) (5.29) 

Assuming =: q and t( 0 ) = 1 , a likely initial value of q( 0 ) is 

0 * 5 TT 5 * Thus the starting values for the nonrotating cantilever beam 
are 



0.00 

0.5773 

0.7952 

0.5773 

-0.7952 


(5.30) 


2 U 

and = 0 , ~ 12,36 . The value u(0) is arbitrarily 

set to zero since it does not matter for a nonrotating beam. Further, 
the control Eq. (5.I8) and the boundary conditions at x = 1 make 


t =0 at x=l ,a fact causing numerical problems , This is 
handled by introducing an end constraint on depth. For the con- 


straint t2:t. (O^x^l) 

min ' ' 


logic is introduced into the pro- 


gram such that in the calculation process whenever t is less than 
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or equal to t , , then t(x) is set equal to t . and St is 

set to zero. The value of St must be zero because no variation of 
depth is permitted, Again^ the integration process continues with 
these new values of t and St . Initially, a large value of t . 
is used, and subsequently it is reduced step-by-step until it is 
*’min ~ • This algorithm is designed to reduce e by one-half 

every time until the depth distribution is zero somewhere between 
X = 0 and x = 1 , and starts again from the preceding step with a 

little higher value of this e . After establishing the likely value 
of £ , convergence is obtained within h to 6 iterations . The dif- 

ferential equations are integrated using the Runge-Kutta numerical 
scheme. The actual depth distribution for the nonrotating cantilever 
beam is shorn in Fig, 5 . 1 . 

Next, a small amount of rotation is introduced, and the corres- 
ponding is estimated from Fig. U.l. Keeping the two param- 
eters and fixed, and using the stationary cantilever re- 

sults corresponding to the optimum depth distribution as starting data, 
the iterative process is repeated to determine the optimum depth dis- 
tribution. This process is continued with a step-by-step increase in 
Vj^ . Figure 5.2 shows the optimum depth distributions with as 

a parameter, 

5.6 PHYSICAL INTERPRETATION OF THE OPTIMALITY CONDITION 

Even in the presence of rotation Eq. {5.2) is self-adjoint, and 
hence the relation [Eq. (5.28)] between the state variables w , p , 
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q , r and the adjoint variables \ \ \ \ is valid • In 

w p q^r 

fact^ it can be easily verified from Eqs, (5*^)^ (5*9)^ and 

(5«19)» However^ this relation is not possible if we include the 
equations 


u"" = - tx 

K = - \vlc\/t^ ( 5 . 31 ) 

in the above sets [Eqs . ( 5 . 26 ) and (5.27)1 in addition to the terms 
due to rotation. Nevertheless, the equality [Eq, ( 5 . 28 )] is valid. 
Therefore, we can write 

^ 5 . 52 ) 

where is a constant* Now^ combining Eqs, (5.18) and (5*52)^ the 

optimality condition can be simplified as 

3 w t - w (Vj^)j^c = Cj, (1 - - c v^ww'x ( 5 . 35 ) 

Integrating the second equation of Eqs. (5 .31), we obtain 

2 h X 

^ f ^2 h ,, 2 '^1^ ^ r ,2, 

*/*^i — 2 “ / '' (5.3^) 

‘^l 0 
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Substitution of Eq. (5.5^) into Eq. (5-55) leads to 


X 

^ ,,2 2 2, v2 ^4 2 2k r , 2 , 

5w t - w X f w dx (5*35) 

0 

Multiplying by t and integrating between 0 and 1 , we obtain 

1 1 

"2, , .2 4 r ^ 2^ 

5t w dx - 

0 0 

11 1 
+ f i f |td| j w^'^dx = f c?tdx (5*36) 

0 ' X / 0 

where ^ is a dummy variable. One can easily recognize that the first 
term is peak strain energy, the second term is peak kinetic energy, and 
the last term is peak potential energy of the equilibrium centrifugal 
forces . 

Thus, the optimality condition [Eq. (5.55)] means that the optimal 
structure is characterized by the property that a linear combination of 
the peak kinetic energy density, the peak strain energy density, and the 
peak potential energy density of the equilibrium centrifugal forces is 
proportional to the design or optimization variable. 

5.7 DISCUSSION OF RESULTS 

The mass saving for a nonfotating cantilever beam is 0.67 at a 
minimum depth constraint of 0.01; on extrapolating for zero minimum 
constraint, the mass saving is 0.68. This number appears repeatedly 



in these sorts of optimizations; indeed, it is interesting to compare 
it with the results reported for a sandwich beam with 50 ^ structural 
mass in Fig, U.lj- of Ref, 5*9* Th® mass saving for this beam is 0.3^. 

A linear extrapolation of this for a solid beam (100^ structural mass) 
gives a mass saving of 0.66, which surprisingly coincides with the 
presently obtained value. 

An increase in rotational speed increases the slope of the depth 
distribution as well as mass savings. Hence, quite sizable mass sav- 
ings can be realized when only one vibration frequency is held constant. 
This is of particular engineering interest since this type of beam is 
used in some helicopter rotor configurations. Thus, if there are no 
constraints such as strength, aerodynamic shape, and blade flutter 
speed, great savings can be realized. 


- 105 - 



VI. THE EFFECT OF ROTATION IN AEROELASTICITY : THE 


WHIRL FLUTTER PROBLEM IN PROP -ROTORS 


6.1 INTRODUCTION 

■flie present chapter extends the study of the effects of rotation 
to aeroelastic problems. Virtually every VTOL aircraft configuration^ 
and its base configurations (the helicopter^ the propeller aircraft^ 
and the turbo-prop aircraft) make use of propellers or rotors for 
thrust or lift. All these configurations present a rich source of 
problems involving the effects of rotation. One potential problem 
that demands careful attention is prop-rotor whirl flutter, which can 
occur in a flexibly supported aircraft -engine -prop -rot or combination. 
It involves coupling of elastic. Inertial, gyroscopic, centrifugal, 
damping, and aerodynamic forces . 

First the vrtiirl flutter problem is briefly reviewed, and then, 
guided by the experience of the studies in the earlier chapters, the 
effects of two additional parameters - steady-state coning angle of 
the blades, and internal damping at the blade flapping hinges — are 
considered. The engine-nacelle-prop-rotor combination is represented 
by a concentrated mass and stiffness mathematical model. Although the 
continuum mathematical model is appropriate for the blades, the con- 
centrated mass model is used in order to gain quick engineering 
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iiisight into this complex problem without extensive computation. With 
this model, all three t 3 T>es of rotors - the first, in t^ich the blades 
are attached to the hub by hinges, free to flap up and down; the second, 
in which the blades are rigidly interconnected to a hub, but with the 
hub free to tilt with respect to the shaft; and the third, in which the 
blades are connected rigidly to the shaft — can be considered. The 
hinges of the freely flapping rotor may be located at varying distances 
from the axis of rotation. 

governing equations of motion with these two additional param- 
eters included are derived in Appendix A by the Lagrangian approach 
using quasi-steady blade element theory for aerodynamic forces. The 
linearized analysis of these equations is applied to two wind tunnel 
models of previous investigations, and the effects of these new param- 
eters on whirl flutter stability boundary are studied. The present re- 
sults are compared with the earlier ones. It is shown that internal 
damping at the flapping hinges may help to resolve the disagreement 
between theory and experiment in some earlier investigations, and to 
explain why some investigations show good correlation between theory 
and experiment and some do not, 

6.2 DISCUSSION OF THE PREVIOUS WORK 

Reed [6.1] covered the early history, the mechanism, and the 
state-of-the-art of the vAiirl flutter problem and showed the exis- 
tence of some gaps, particularly in correlating theory and experiment 
in the case of a flapped blade prop-rotor in the propeller mode of 
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operation. Comparing analytical and experimental results for two flap- 
ping hinge offset cases, one 15 ^ and the other 8^, with Isotropic pylon 
stiffnesses, he found that theory and experiment are in good agreement 
for the hinge offset case. However, the experimental model for the 
8'^ hinge offset case fluttered in a forward whirl mode, and this forward 
whirl flutter was not predicted by his analysis in spite of introducing 
aerodynamic lag effects as high as 1 + 5 °. Earlier to this work, a simi- 
lar type of disagreement between theory and experiment was reported by 
Richardson and Naylor [ 6 . 2 ]. This work was followed by an article by 
Young and LytvTyn [6.5], who studied analytically the influence of flap- 
ping restraint on the whirl flutter of a fully flapped rotor (a rotor 
with only flapping hinges without any offset). They established the 

existence of an optimum flapping frequency parameter between 

pi o 

1.1 and 1.2, depending on the ratio of the nacelle (pylon) stiffnesses 
in the pitch and yaw directions, and that at the optimum condition 
nacelle stiffness requirements for neutral stability are at a minimum. 
Further, it was shown that forward whirl type instability is possible 
if the pylon mounting is anisotropic, or the flapping frequency param- 
eter is less than the optimum value. 

Almost at the same time as Ref. 6,5, motivated by experience in 
conducting an exploratory investigation of whirl flutter in high by- 
pass ratio fan jet engines (Refs. 6 . 1 + and 6.5), Kaza [6.6] studied 
both analytically and experimentally the whirl flutter problem in 
prop-rotors for two hinge offset positions ( 10 ^ and 13 . 6 ^). Both 
theory and experiment predicted backward whirl type instability for 
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both the hinge offset cases with isotropic pylon mounting, and the 
quantitative agreement is reasonably good. Also it was shown that an 
optimum flapping frequency parameter, similar to that of Ref. 6 . 5 , ex- 
ists through a combination of hinge offset and flapping spring restraint 
parameters. The possibility of both backward and forward whirl flutter 
is established when the pylon mounting is isotropic and the flapping 
frequency parameter is less than approximately I. 05 . Moreover, the 
dominant mode is backward whirl in the range 1.05 to 0 . 95 » forward 

whirl in the range 0.95 and less, of the flapping frequency parameter 
for isotropic pylon mounting. Thus, this study indicates the possibil- 
ity of forward whirl even when the system is mounted isotropically. 
However, both this study and that of Ref. 6,3 are unable to defend the 
experimental observations of Ref. 6,1 for the 8 ^ hinge offset case and 
those of Ref, 6,2. 

In attenq>ting to explain the disagreement between theory and ex- 
periment for the 8^ hinge offset case, Gaffey [6.7] assumed that the 
flapping hinge also acts as a lead-lag hinge. This assumption opened 
the room for the possibility of forward whirl instability for both the 
8^ and the 15 ^ hinge offset cases, and surprisingly the theory and ex- 
periment for the Qio case compare well . However, the forward whirl in- 
stability due to the lead-lag motion for the 13 ^ case was not at all 
encountered in the wind tunnel. Evidently, the disagreement was still 
unresolved. Recently, Baird, Bauer, and Kohn [6.8] presented results 
of a parametric study on whirl flutter. This study includes vari- 
ations of blade pitch-flap coupling, flap hinge offset, pylon damping. 
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and pylon stiffness. A good agreement between theory and experiment was 
reported for both the 5^ and the hinge offset cases , 

More recently, almost at the same time as the present work, Kvaternik 
[6.9J presented results of experimental studies on whirl flutter, those 
of joint NASA/Sell Helicopter investigations of a 0 .IJJ-scale, semispan, 
dynamic, aeroelastic model in the NASA-Langley Transonic Dynamic Tun- 
nel, and a comparison with theoretical results. Analytical results 
were shown to be in good agreement with measured dynamic character- 
istics , Kvaternik also covered the mechanism of whirl flutter in prop- 
rotors and gave a comprehensive list of references on the subject. His 
study also included a small steady-state coning angle of the blades arid 
a small damping at the flapping hinges . Confining his model to within 
the covered range of these parameters, Kvaternik reported that an in- 
crease in coning angle has a destabilizing effect, and an increase in 
flap damping has a stabilizing effect on whirl flutter. Though this 
study shows good correlation between theory and experiment, apparently 
it does not aid in resolving the disagreement between theory and experi- 
ment noticed in Refs. 6.1 and 6 . 2 . 

6.3 STATEMENT OF THE PRESENT PROBLEM 

The problems considered in the previous chapters suggest that damp- 
ing may be a factor for the disagreement between theory and experiment. 
Also, the results of Ref. 6.6 show the possibility for predominant for- 
ward whirl flutter if the flapping frequency parameter is less than ap- 
proximately 0,95, which is possible according to Section 3.7 if there 
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exists a steady-state coning angle of the blades. In fact, a small 
coning angle results in practice from the balance between centrifugal, 
gravity and aerodynamic forces at the unperturbed steady-state condi- 
tion of the rotor in the case of flapped blade systems, or from the 
steady flap bending in the case of rotors with flexible blades attached 
rigidly to the hub. Therefore, in the present study this parameter is 
also considered. Also, Loewy [6.10] pointed out that the coning angle 
may be the source of error for the discrepancy between theory and ex- 
periment. 

The problem of this chapter is to study in depth the influence of 
these new parameters on the whirl flutter boundary. The contents of 
this chapter are summarized in a paper by Kaza [6.11], 

6,h MATHEMATICAL MODEL 

The discretised mathematical model shown in Fig. 6,1 consists of 
N (N s 5) ideally rigid blades, restrained elastically with equal 
spring constants at the flapping hinges, A viscous type of damping, 
which may be due to friction at the flapping hinge in the case of 
flapped rigid blade rotors, or may be due to hysteresis type damping 
of the blade material in the case of rotors with flexible blades rig- 
idly attached to the hub, is considered. This rotor system along with 
an ideally rigid pylon or nacelle is restrained elastically in pitch 
and yaw at a distance h from the rotor plane. A viscous type of 
damping is also assumed in the pitch and yaw directions. This dy- 
namic system has N + 2 degrees of freedom with respect to a rotating 
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observer. The rotor blades also have a steady-state coning angle , 
which is due to the balance between centrifugal, gravity, and aerody- 
namic forces at the unperturbed steady-state position of the rotor in 
the case of flapped blade systems or due to a built-in angle, in the 
case of rotors with flexible blades attached rigidly to the hub. If 
N ^ 3 j the rotor has rotational symmetry. Consequently, the N de- 
grees of freedom of the rotor blades may conveniently be reduced to two 
degrees of freedom through quasi-coordinates, namely, the pitch flap 

angle p and yaw flap angle p . Since attention is directed to 

^ \|r 

the propeller mode of operation of the prop-rotor, the resulting equa- 
tions of motion are free from time dependent coefficients. 

Quasi-steady blade element theory is used to generate the aerody- 
namic coefficients . A typical blade element and velocities experienced 
by it are shown in Fig. 6.2. The quasi-steady blade element theory is 
justified mainly because the vorticity shed by the blades is carried 
away from the propeller plane at a much faster rate than that associ- 
ated with the helicopter in forward flight. In addition to wake ef- 
fects, several other assumptions, which are quite common in whirl 
flutter stability analysis, are made. They are: 1) no compressibil- 

ity effects; 2) no aerodynamic forces on pylon and spinner; 3) the 
blade airfoil is symmetric with constant chord along the span; 4) 
lift coefficient is 27T and profile drag coefficient is a quadratic 
function of angle of attack; and 5) induced velocity due to thrust. 
All these assumptions are self-explanatory except the last one, which 
is justified since the disc loading in the propeller mode of airplane 
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mode of operation is less than that in the hover mode. No restriction 
is imposed on the value of the steady-state coning angle 

Finally, the governing equations of motion of this mathematical 
model are derived in Appendix A. They are a set of four simultaneous, 
ordinary, coupled second order differential equations with respect to 
a space fixed coordinate system. 


6.5 MATHEMATICAL ANALYSIS 


The system with the blade motion & , and P fixed, and with 

the coning angle set equal to zero, is considered. This has two 

degrees of freedom, pitch and yaw. These are uncoupled in the absence 

of rotation. When the rotation is introduced, these two motions are 

coupled. The higher frequency one is called forward whirl, and the 

lower frequency one backward whirl . 

When the pitch and yaw freedoms of the pylon are locked, the rotor 

blade has an undamped flapping frequency CD , and a damped frequency 
— , P 

1 - L with respect to the rotating observer. In the 


Pld ■ pi 


absence of steady coning angle and damping , the flapping 

frequency varies from to infinity. The lower limit corres- 

ponds to an ideally rigid blade without an elastic restraint at the 
flapping hinge, which is located on the axis of rotation, and the upper 
limit corresponds to an ideally rigid blade attached rigidly to the hub. 
But in the presence of the coning angle p^ and the damping , the 

lower limit of will be less than , depending on the value 

of Pq • When the blade motions are expressed with respect to a space 
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fixed coordinate system, the flapping frequency is split into 

two frequencies, ■ ^ ■ “pid and ^ frequencies occur 

in pairs as expected from a gyroscopic system. The rotor tip path 
plane has two types of motions, depending on the parameters, with re- 
spect to a stationary observer. When ^pXd ’ positive 

frequencies are - H and . The first one corres- 

ponds to a backward whirl mode and the second one to a forward whirl 
mode of the rotor tip path plane with respect to a stationary observer. 

For the case > cu , , , the two positive frequencies are H , 

o pld ^ '1 o pXd 

and , both representing the forward whirl modes of the tip 

path plane of the rotor. Thus the presence of the coning angle and 
damping may even change the rotor model characteristics, depending 

on the other parameters . 

The governing Eqs . (A,3T) of motion of the model in a uniform, in- 
compressible, and inviscid flow at a speed V are derived in Appendix 

A, 

Assuming 


[q] 



+ i7)T 
e 



( 6 , 1 ) 


Eqs. (a. 37) can be written as 


[D] (Q) . . i CQ} 


( 6 , 2 ) 


- 115 " 



where 


[D] = 

tv 

[D^] 


-[1] 

[o]_ 




[D^] = 

[[SS] - 

[AS]]'^ 

[D^] = 

[[SS] - 

(q. ) 

[AS]]'^ 

[Q] = ' 

1 o J 



[AN] 

[SD] - [AD] 


( 6 . 3 ) 


This is a standard eigenvalue problem. The eigenvalues are four com- 
plex conjugate pairs, two pairs for the pylon and two pairs for the 
rotor. The general solution of the system will be obtained if we con- 
sider the eigenvalues only with positive imaginary parts. The associ- 
ated eigenvectors will determine the direction of the whirl mode. When 
the flapping frequency o£ the blade with damping and aerody- 
namics included is greater than , there are two backward whirl 

modes with frequencies m = (pylon backward) and o>_ = 

1 1 o ' 3 5 o 

(rotor backward), and two forward whirl modes with frequencies 
"^2 = ^2% (pylon forward) and (rotor forward). But when 

the backward rotor mode changes to forward whirl, and there 
will be three forward whirl modes and one backward whirl mode. The 
system is critical when the real part (p^) of any one of the eigen- 
values is zero. Keeping all the parameters, e , p j Co > ^ ^ f 

o p Q ^ 

^ '^h * ^ and c as constants, is varied 
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until the system is unstable in one or possibly more than one mode. The 
stability boundaries are obtained using the complex eigenvalue algorithm 
(Ref. 6.12). 

As mentioned earlier, this analysis is applied to two models, the 
and 13 ^ model of Ref. 6.1, and a fully flapped (e = o) version of 
the model of Ref. 6.6. Most of the parameters of these models are taken 
from the appropriate references. Some parameters for the 8^ model of 
Ref, 6,1 were obtained through private communication from Reed, The 
additional parameters needed to consider the effect of coning angle 
are calculated from the available data. Finally, all these parameters 
are tabulated in Table 6.1. 


6.6 RESULTS AND DISCUSSION 


The main attention is focused on the effect of and on an 

isometrically mounted pylon system, which is the most critical case from 
the flutter point of view, is also included since the effect of 

and Cp depends on it. 

The effect of coning angle and flapping restraint stiffness on 
flutter speed is shown in Fig. 6.3 for the model of Ref. 6.6. Depend- 
ing on the system parameters, either whirl flutter or static divergence 
may be possible in all except one of the modes, the rotor high frequency 
forward mode. However, which mode is predominant depends on the param- 
eters , The magnitude of the critical speed strongly de- 
pends on both and , For each coning angle , there 

exists an optimum flapping frequency ratio at which the flutter 

p y 
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TABLE 6.1 


Farametet's of the Models Used in Whirl Flutter Calculations 
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TABLE 6.1 (Cont*d) 


Parameter 

Fully Flapped 
Model, 

Ref. 6,6 

8^ Case 
Ref. 6.1 

15 fa Case 
Ref. 6.1 

a 

0 

0.0 

0.0 

0.0 


varied 

varied 

varied 


varied 

varied 

varied 


0 .0087 

O.OO8T 

0.0087 

1—1 

-0 .0216 

-O.02U6 

-0.0216 

^iy2 

o.k> 

O.I4O 

O.iK) 
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speed is maxiinum, and below that value there is a possibility of forward 
type whirl flutter even when the system pylon mounting is isotropic. At 
the optimum condition, the flapping frequency of the rotor with respect 

. p , p jL 

to a rotating coordinate system is to = ft [cos 2p + (to /to ) (to /ft ) ]^, 
where is the associated flutter speed. In the present case, this 

value is approximately 1.05 ^ for all values of P considered, and it 
is in agreement with those of Refs, 6*5 6 , 6 . Hence^ it appears that 

the optimum condition can be achieved either by changing p or 

o p 9 

However, an increase in P causes a decrease in flutter speed \dien co 


61 


is less than 1.05 p increase in flutter speed if m is 

o pi 

greater than 1 .05 . In view of the fact that both the flutter speed 

and its mode are very sensitive to in the range 1.00 ^ to 1.1 

pi o o 

which, in turn, depends on p , one must include p in the analysis* 

even for the windmilling case, where is minimum. In fact, the 

models tested in Refs. 6.1, 6.2, and 6.6 have in the range 1,05ft 

pi o 

to 1.20 ft^ without any effect. Figure 6,k is a cross plot of Pig. 

6,3, showing explicitly the effect of on both flutter speed and 

mode for the case of a fully flapped rotor without an elastic restraint 


spring at the flapping hinge and with isotropic pylon mounting. Both 
pylon backward and forward modes flutter. The backward mode is pre- 
dominant between 0*^ and 15 ° of coning and the forward mode takes the 
lead when the coning is more than 15°. The overall effect is a de- 
crease in flutter speed with an increase in coning. 

Next, the effect of coning for the case of a hinge offset blade 
rotor is considered. The parameters used are those of the 8^ hinge 
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offset model with windmilling blades of Ref. 6,1. The range of is 

extended to the negative side in order to account for the windmilling 

case. The flapping frequency of the rotor with respect to a rotating 

o o JL 

observer is given by m -ft [cos + cos p 1 /ir. + 

Figure 6.5 shows the effect of coning angle on the aforesaid model with 
out elastic restraints at the blade hinges . An increase in from 

-5° to about causes an increase in flutter speed, and from there on- 
wards it not only causes a decrease in flutter speed but also destabil- 
izes the pylon forward mode in addition to the backward mode. However, 
the predominant flutter mode is backward except in the range 6 ° to 8 °, 
which corresponds to = 1.05 . Studying closely Figs. 6 ,k and 

6.5, one can quickly conclude that the effects of 3^ , e , and 

can be integrated as the effect of co , if p is positive. But when 

pi o 

the propeller is windmilling, is negative, and hence its effect 

should be considered explicitly. At the optimum condition both diver- 
gence and flutter occur almost at the same speed, as in the case of a 
fully flapped rotor. This model exhibited a forward whirl instability 
in the wind tunnel. Both the analysis of Ref, 6.1 and the present one 
for = 0 predicted a backward whirl flutter, and they are in good 
agreement. This model in the windmilling case might have a coning 
angle anywhere in between -5° and - 2 ^, In this range there is no 
possibility for forward flutter as expected from experiment. However, 
a forward whirl flutter at a much higher speed than observed in the 
wind tunnel is possible when is in the vicinity of 7 . 5 °. On 

further increasing to 15°, the predominant flutter mode changed 
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from forward to backward. The magnitude of this backward speed is in 
agreement with that of the forward speed noticed in the tunnel. Thus 
the effect of the coning angle may not be the reason for the disagree- 
ment between theory and experiment on this models but it certainly has 
a marked influence on the stability, both divergence and flutter. 

The effect of on the flutter speed for the case of fully 

flapped rotors with as a parameter is shown in Figs. 6.6(a) 

and 6.6(b) for two values of . The damped flapping frequency of 

the rotor is given by ^ 1 - , with respect to the blade 

fixed axis system. It is observed that a small amount of damping ^ 

(^^ of critical) at the flapping hinge not only changes the predominant 
flutter mode from backward to forward but also reduces the flutter speed 
This effect is analogous to that of internal damping in the rotor as dls 
cussed in Ref. 5.9. The difference in flutter speed between the pylon 
backward and forward modes decreases with an increase in the elastic re- 
straint parameter and finally changes the predominant mode to 

backward whirl as expected. 

Extending this study on the effect of flap damping to hinge offset 

bladed rotors, the parameters of the 8^ and the hinge offset models 

are selected. Unfortunately, the data on damping at the flapping hinge 

are not available, and hence the flutter boundary ^ versus H /ou 

P o Q ’ 

with as a parameter, is shown in Fig. 6.T. Pylon damping effects 

in the presence of flap damping are also included in Fig. 6.7(b). As 
in the case of fully flapped rotors, the damping at the flapping hinge 
continuously reduces the backward flutter speed. It might even 
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destabilize the forward whirl mode at a much lower speed than that of 

the backward whirl mode when is more than 0.l6. Also^ when Co = 0 , 

P P 

pylon damping has no appreciable effect on flutter speed. When is 

I 

more than 0.12^ pylon damping has a very little effect on backward whirl 
mode instability and has a marked stabilizing effect on forward v^irl 
instability. Comparing the present value of the flutter speed, /to 

O C7 

= 8.^ at Cq=0,P =0 , with the corresponding one (fi /to = 8.^4^) 

p O O C/ 

evaluated from Fig. 9 of Ref. 6.1, shows that they are in good agree- 
ment. However, the present theory does not account for aerodynamic lag 

effects whereas Ref. 6.1 does. But when Co ^ 0.1 , and C^ C 

P y \[f 

lie in the range 0 to 0.04, the predominant flutter mode is forward and 
the speed Q /to lies in the range 1.7 to 5*5* Hence, the forward 
flutter instability observed in the wind tunnel within the range 0 /cjci 

O u 

= 1.52 and 2,^4- (from Fig. Ref. 6.1) may be due to damping at the 
flapping hinge. 

Figures 6.8 and 6.9 summarize the calculations for the IJ^ hinge 
offset model. The variations of both real and imaginary parts of the 
eigenvalues with 0 /to are shown in Fig. 6.8 along with the corres- 
ponding ones of Ref. 6.1. Comparing these results, the agreement is 
very good. The difference in the variation of damping for the rotor 
backward whirl mode may be due to the fact that we assumed an aerody- 
namic phase lag of zero whereas a lag of 15 ° was assumed in Ref. 6.1. 

The flutter mode is backward rotor ^irl as observed in the wind tunnel . 

The next question concerns the effect of C^ on the flutter speed 
for this case. The answer is in Fig. 6 . 9 , which shows the variation of 
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FIG. 6. 8- -Calculated flapping bla( 
0> /cD =1.0 
6 ®- 


0) = 0 * f 


1 PYLON BACKWARD 

2 PYLON FORWARD 
5 ROTOR BACKWARD 
h ROTOR FORWARD 




e whirl modes — e = O.IJT 
= 0 ; ^0=0 (see Table 



as a parameter. It can be seen 


flutter speed with , with 

P y \|f 

that there is no possibility for forward whirl flutter for in the 

range 0 to O.52, unlike the 8^ hinge offset case. This reversal in the 

^ effect is due to an increase in the flapping frequency parameter 
P 

CO-. . which increases with the hinge offset parameter e . This is 
pi 

precisely what is expected from the calculations of a fully flapped 
rotor for large values of the flapping restraint parameter > 

see Fig. 6.6(b). From Fig. 6.9, again it is clear that pylon damping 
has no appreciable effect on flutter speed. 

In summary, it is seen that internal damping and coning angle have 
a marked influence on whirl flutter stability, depending on other param 
eters. In fact, the internal damping at the flapping hinge not only 
changes the flutter speed but also its predominant mode, depending on 
the other parameters. It is believed that a proper account of this 
internal damping might help to explain why some investigations show 
good correlation between theory and experiment and some do not. 
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VII. CONCLUSIONS 


This dissertation covers four areas, vibration, stability, opti- 
mization, and aeroelasticity, considering elastic continuum models and 
concentrated mass and stiffness models. This study is undertaken with 
the following objectives: a) to group together at least some, if not 

all, of the problems involving effects of rotation; b) to identify the 
common features of this class of problems; c) to use this experience 
and observations profitably to study some problems which have not yet 
received a special treatment of their own, such as the stability of a 
coned beam and the least weight design of a rotating beam with a con- 
straint on fundamental flapping frequency; and d) extrapolate this in- 
formation for better understanding of the xdairl flutter problem in 
flapped blade prop-rotors, particularly to resolve existing disagree- 
ment between theory and experiment in some earlier investigations under 
certain conditions. 

At the outset, it is realized that the continuum mathematical 
model in the presence of rotation has a practical limit to the class 
of problems for which analytical solutions are possible. Adapting the 
common engineering philosophy that the understanding of any problem is 
helped if its relation to other problems is appreciated, the equations 
of motion of a general unrestrained, deformable, and rapidly rotating 
body are developed. The explicit need for the nonlinear theory of 
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elasticity under certain conditions is emphasized. Some common fea- 
tures such as the choice of axis system^ self-adjointness^ the phenom- 
enon of frequency splitting, shortcomings of stability methods as ap- 
plied to gyroscopic systems, and the consequent need for introduction 
of imperfections to extract all the instabilities by kinetic method of 
stability, are identified and discussed. 

On the basis of a linear theory, vibration and stability of a ro- 
tating coned beam is examined. It is seen that the beam has several 
interesting bands of instability depending on the coning angle and the 
rotational speed. It must be pointed out that we assumed that the 
bearings are rigid and axially symmetric. If the beam is unsymmetric, 
and the bearings are considered to be flexible and also unsymmetric, 
the equations of motion will contain periodic coefficients. The solu- 
tion of such partial differential equations is not at all possible at 
present. Numerically exact frequencies of the beam are obtained, and 
these results are used to validate approximate methods. 

The field of parameter optimization and optimal control theory 
techniques are briefly reviewed, and a start is made in applying these 
techniques to simple rotating systems. Though these continuum models 
have practical limits as stated earlier, it is believed that these re- 
sults can be extrapolated to understand the behavior of complicated 
systems. Optimum depth distributions of a solid rotating cantilever 
beam are obtained with rotational speed as a parameter and with a fixed 
first flapping frequency. It is shown that, even in the presence of 
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rotation, the optimal structure is characterized by the condition that 
Lagrangian density is constant. 

Finally, guided by the experience of the above general study, the 
effect of two parameters — steady coning angle of the blades and in- 
ternal damping at the flapping hinges — on whirl flutter in prop-rotors 
are studied. It is shown that these two parameters have a marked in- 
fluence on stability. The effect of this internal damping for a cer- 
tain range of model parameters is discovered to be analogous to that 
of internal damping in gyroscopic systems. It is also shown that this 
internal damping may help to resolve the disagreement between theory and 
experiment in the earlier investigations under certain conditions. 

The study on the effect of rotation has not been completed in its 
entirety but several useful observations have been made and the exper- 
ience gained has been applied to some practical cases. There is a vast 
amount of work to be done. Some examples are the dynamic response of a 
rotating coned blade with random inputs, the influence in depth of in- 
ternal and external damping in gyroscopic systems, and continuation of 
optimization studies in rotating systems with aeroelastic constraints. 
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APPENDIX A 


EQUATIONS OF MOTION USED IN WHIRL FLUTTER STUDY 


The mathematical model of an idealized dynamical system considered 
here is shown in Fig. 6.1. The governing equations of motion are de- 
rived with respect to an inertial coordinate system using the classi- 
cal Lagrangian approach and quasi -steady blade element aerodynamic 
theory. The coordinates of a point on the blade with respect to the 

inertial axis system X , Y > Z ,0 can be written as 

s s s s 



(A.l) 


where 




COS 9 0 

0 1 

-sin 9 0 


sin 9 
0 

cos 9 


(A.2) 
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and 



cos ^ -sin i|f 0 

sin !• cos ijf 0 

0 0 1 



~1 
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0 

[T*l - 

0 

cos 

d> 

n 

-sin 

n 


0 

sin 

4* 

n 

cos ^ 

n _ 


cos 


0 

sin P 

n 


0 


1 

0 


1 -sin 

0 

cos P 

n 


Aligning the axis along the axis of the blade, 
position vector of an arbitrary point on the axis with respect 
blade fixed axis system is 



(A .5) 


(A.M 


(A. 5) 

the 
to a 


(A. 6) 
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When the hinge offset of the blades is uniform, the position vec- 

tb. 

tor of the flapping hinge of the N blade is 



(A.7) 


Considering the case that the pylon is restrained elastically in pitch 
and yaw directions at an equal distance from the propeller center, we 
can write the position vector of the point 0 as 



(A.8) 


Analogous to , the position vector of an arbitrary point on the 
axis of the symmetric nacelle with respect to the inertial axis sys- 
tem becomes 



(A. 9) 
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The expression for the total kinetic energy of the system is 


N 


R 


T = 


-. 1 ./ 


tp •- 


dm^ 


1 / 


{x 3^ {x 3 

gn-’ '• gn-’ 


dm 


n 


(A. 10) 


n=l 


0 


0 


All the ingredients required in this equation are given by Eqs . A.l 

to A. 9. The flapping angle of the blade is expressed with 

respect to the rotating coordinate system and 0 

However, this can also be expressed conveniently for the case of rigid 

blades with respect to the inertial coordinate system X , Y , Z 

s s s 

and 0 through quasi -coordinates and p and a steady-state 

s t/ \|f 

coning angle as 




n 


P. + P. 


n 


Po + cos 


•J' +6 sin <t> 

n . n 


(A.ll) 


This simplification is possible because the other possible patterns of 
the blade motion do not couple with the pylon freedoms. Considering 
6 as small* we can further write 

cos P 

cos p = cos p - p sin p - p 

n o n o n 


_ sin P 

sin p = sin p + p cos p - ^ p 

n o n o ^ n 


(A. 12) 


Since 9 and ^ are small^ the trigonometric functions of these co- 
ordinates can be expanded in power series, and the terms higher than 
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the second order can be neglected. The azimuth angle for N 

equispaced blades can be written as 

% = <» + I (N-l) (A. 13) 

where 

<l> = Ht (A.l^) 

Substituting Eqs . A, 11, A, 12, A. 15 and A. 1^1^ into Eq, A. 10 and 
applying Lagrange’s equation (Ref. 2.1), we obtain the following gen- 
eralized inertial forces: 

Pg = {V" 1 ^; * 

+ [Ip,(l + cos 2p ) + 21 cos 6 ] p'} (A. 15 ) 

P o o o ^ 

- [Ip(l + cos 2P^) + 21^ cos P^] P'} (A. 16 ) 

p. = {I p! + 21 p' + I e" 

Pq P 9 p i|r 2 

+ [Ip (cos 2p^ - 1) + I^ cos p^] p^] 

+ [Ip(l + cos 2p^) + 21^ cos p^] p' (A. IT) 
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and 


where 


- IVl * hi ■ ^V'e 




+ p [(l„ cos 23 - 1) + I cos 3 ] 


- [ 1^(1 + cos 23 ) + 21 cos 3 ] e'} 

3 O O O 


(A. 18) 
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dr 
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J* m^r r'dr 
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s ds 
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Iq + I COS 3 + — ^ sin 3 

P o o 2 ° 




= Nh 
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dr 


(Cont’d) 
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and 


M = N ^ 
0 


Mh^ = Mh^ + 21, sin g + 2I„ sin^ P 
no o (3 o 


1^ = Il^Mh 


I- = (I. - 21 ) + 21 cos p 

0 ^6 o' o o 


+ I„ sin"^ p + 2L sin 6 
6 no 


I = I, + Mh + I„ 
\|f 1 N 


I = (I - 21 ) + 21 cos p 
f llf o 


+ Ip sin P^ + 21^^ sin p^ 


T = ft t 

O 


(y - r n 


(Cont’dj 
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and 



(A. 19 ) 


yaw, 


The potential energy due to linear elastic springs in the pitch, 
and flapping inodes can be written as 


N 


U. 


PE 


i / + 


2 


r 






(A. 20 ) 


n=l 


Substituting Eqs . A. 11, A. 12, A. 13 and A.l4 into the above equation, 
we get 




(A. 21) 


Then the generalized elastic forces follow from Eq. A. 21 




(Cont’d) 
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and 


where 


f 


T 2„2 
17^ \|f 
^ 'I' 




2 2 


0 


i 


2 2 

V^‘' % 


(A. 22) 
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fl 


CD 

_g. 

ft 



ft 




(A. 23 ) 


A viscous type of damping is assumed in the Q , f , and p 
motions . Then the dissipation potential is given by 
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N 


1 ‘2 1 *2 1 ‘ 


*2 
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(A.2U) 


n=l 
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Again the Lagrangian equations lead to the following generalized damp- 

i 

ing forces 


where 






7 i V 

o i(r \)f 








21 05 
t|f i 








^c' 
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'pi 


f I \ 2 

cos 2p + — cos P + I 
o T 'p ; 

^ ip / 


(A.25) 


(A.26) 
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The quasi -steady aerodynamic forces are derived from the well 
known blade element theory of the propeller* A typical blade element 
is illustrated in Fig. 6.2. The elemental lift and drag expressions 


are given by 


dL^ = I pcU^C^dr 

dD = ^ pcU^C ,dr 

n ^ d 


where 




= K 


do 






By resolving dL^ and dD^ 
the elemental inplane force dH^ 


, the elemental thrust dT 

n 

can be obtained as 


dT 


n 



- UU 


K 


do 


+ K 


dl 


^V-tan 



(A. 27) 


(A. 28) 
and 
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1 

— pcdr 

2 



/ -1 
I V-tan 


U \ 

"t / 


+ mj„ 


^do " ’'dl 


v-tan 



(A.29) 

where U and U_ are the perpendicular and inplane components of 
P T 

velocity in blade fixed coordinates* 

Referring to Fig. 6*1^ the velocity vector with respect to the 
inertial axis system is given by 




(A .30) 


From Eqs . A,1 and A .50 the velocity components experienced by a 
blade element can be written as 





{V.) (A. 51 ) 
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Expansion of Eq* A.^1 yields 


U 

P 


+ 


V(cos 6 cos cos - cos 6 sin ^ sin. (p^ sin 

sin 9 cos cp sin p ) + 0 cos ^ cos cp (r + h sin p ) 
n n T n n 


^ sin cp^(r + h sin p^) + r'p^ 


Uj = V(cos 9 sin cos <p^ - sin $ sin 9 ^) 

* 

+ r0(sin ip cos P + cos -jr sin 9 sin P ) 

T n T n n 

» • 

+ h© sin CD cos y\r - cos cp sin p 
n ^ ^ n n 

• • 

- cos cp^ + cos H- r^9 sin ^ (A, 52 ) 

The elemental force vector inertial system follows 

from Eqs. A*1 and A*29 


[dF„] 


(laUT HtJIt.1 

0 i|f 9 P 



(A. 55) 


From the principle of virtual work, the generalized aerodynamic 
forces are 

N R 


■ Z) / 


n=l 0 


„ dx 
T s 

S<Ii 


(A.5^) 
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These forces are nonlinear. Hence, they are expanded in Taylor series 
about the steady initial values of 9 » ) an.d , and the series 

are truncated at the first derivatives , For a particular case in which 
the blades are twisted such that CC^ is constant, the linearized aero- 
dynamic forces, after neglecting the effects of the hub and the air- 
loads on the nonflapping portion of the blades, are given by 


+ MgS' - 

o = + M6^|r} 

t f 

Qe = [(M, - - (M^ -H M' -f 

+ (M^ + - Wp 0 - (Mj^ + + Np ^ 

+ (M^ + P; - (»5 N') p;|^ + (Mg + - N') 6' 

- (M^ + M;^ + Ng + N^) ^'] 


% “ “o '<“2 + «2 + «1 + »I) Pfi + («1 + M' - 

+ (Mj^ + Mf; + + Np 0 + (M^ + M^ - Ni^ - Np ^ 

+ (N^ + Np P' + (M^ + Mp P' + (M^ + + Ng + Np d' 

+ (Mg + M^ - - Np ^'] (A.55) 
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do dO dl o d2 o 


e = r /r 
e' 



= t/r (A.56) 


The final equations of motion are obtained by simply adding to- 
gether all the forces acting on the system, and finally they are 
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written below in matrix notation as 


[AN][q"} + 


[SD] 


[ad] 


{q'} + 


[SS] 


[AS] 


{q} = 


where 


{q} 



[an] 



(A .37) 
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21, 


•P 


(l+cos2P )+2I cosp 
P' o' o o 


0 




Ip (l+cos2P^ )+2I^cosp^ 


21 , 


- 21 . 


•[lp(l+cos2p^)+2I^cosp^] 




■[Ip (l+cos2p^ )+2I^cosP^] 


- 21 , 


2? I 7 

t t i 


Ip(cos2p^-l)+I^cosP^+7pIp 




0 


I 7^ 
e^d 




Ip(cos2^^-l)+IoCos^^-f7pIp 


0 


I 7^ 

t i 


(Cont'd) 





+ M' - (N. 


- Nj^ - Nj; -(M^ + + Np - (Mj^ + + N^ 


+ N^ + N' Mj^ + Mj^ + N^ + N ' 


Ml + Mj - N^ 


M5 + - Nj; 


+ M^ Mg + - N^ - N' -N^ - N' - (M^ + - 1 - Ng 4- N^ 


N^ + N' 


M^ + + Ng + + M' Mg + - N^ - N' 


(A .38) 
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